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(i) 


ABSTRACT 


This  thesis  deals  with  second  order  partial  differential 
equations  of  elliptic  and  parabolic  types  containing  a  small  parameter 
e  associated  with  the  second  order  derivatives.  Consider  problems 
I  and  II  below  with  operators  and  L^,  respectively.  The 

differential  operator  L£  is  given  by  L£  =  eL^  +  where 
is  a  differential  operator  of  the  first  order  and  L  is  a  differential 

i- 

operator  of  the  second  order.  If  is  the  solution  of  problem  I 

having  PDE  L£[u£]  =  f  satisfying  the  given  boundary  or  initial 

conditions  and  u^  is  the  solution  of  problem  II  having  PDE  Lq[uq]  =  f 

satisfying  an  appropriately  modified  initial  condition,  then,  for  e 

small,  we  can  argue  that  the  two  solutions  should  be  quite  similar, 

at  least  in  some  sense.  The  problem  is  to  determine  conditions  under 

which  u  u_  as  e  -*  0.  In  addition  we  examine  how  the  difference 

e  0 

between  u^  and  u^  behaves  on  those  portions  of  the  boundary 
where  u^  does  not  satisfy  the  prescribed  condition  on  u^ .  Indeed 
the  difference  consists  of  two  terms,  the  first  of  which  is  small 
(say  0(en))  uniformly  in  e  in  the  whole  region  of  the  problem 
and  the  second  of  which  is  appreciable  only  near  the  boundary  but 
is  (say)  0(e  )  away  from  the  boundary. 
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CHAPTER  I 


INTRODUCTION 


Suppose  Lq  is  a  first  order  linear  partial  differential 
operator  whose  coefficients  are  functions  of  x  and  y  and 
is  a  second  order  linear  partial  differential  operator  with  constant 
coefficients  of  either  elliptic  or  parabolic  type.  Let 

L  -  eL.  +  Ln  . 
e  10 

Suppose  u(x,y;£)  is  the  solution  of  the  problem 

Vul  ■ £ 

in  an  open  domain  D  subject  to  suitable  boundary  conditions.  If 
we  call  U(x,y)  the  solution  of  the  associated  reduced  equation 

L0[U]  -  f 

over  the  same  domain  D  subject  to  the  same  boundary  conditions  as 
u  on  at  least  part  of  the  boundary,  then  it  has  been  shown  that 
under  certain  conditions  that 

u(x,y;£)  -*  U (x , y )  as  £  -*  0+ 

uniformly  over  most  of  D.  The  regions  in  which  this  condition 
does  not  hold  are  to  be  found  near  the  boundaries  where  the  values 
of  U  and  the  arbitrarily  defined  values  of  u  do  not  coincide. 

To  accomodate  this  difference,  a  boundary  layer  term,  called  z(x,y;£) 


2 


in  this  thesis,  is  introduced. 

There  is  a  rapidly  growing  literature  dealing  with  this 
type  of  problem.  The  general  survey  of  the  work  done  on  this  subject 
before  1955  can  be  found  in  [6]..  However,  in  this  thesis,  we  are 
interested  in  only  the  second  order  elliptic  and  parabolic  equations 
which  can  be  dealt  through  the  use  of  Wasow's  technique  [10]  and 
Levinson's  technique  [7], 

Now  in  the  case  of  the  elliptic  problem,  the  result 
obtained  by  Levinson  [7]  is  of  the  form 

u(x,y;e)  =  U(x,y)  +  z(x,y;e)  +  w(x,y;e) 

where  w  =  0(/e")  uniformly  over  D.  In  this  thesis,  we  show  how  this 
error  or  remainder  term  may  be  refined  to  w  =  0(e);  and  indeed 
present  a  method  of  obtaining  a  representation  of  u  in  D  in 
which  the  remainder  term  in  0(en)  where  n  is  any  positive  integer. 
This  type  of  analysis  is  extended  further  to  the  case  where  the 
associated  reduced  equation 

Lq[U]  =  a(x,y)Ux  +  b(x,y)U^  +  c(x,y)U  =  f(x,y) 

is  replaced  by  a  similar  first  order  operator  in  which  a,  b,  c,  f 
as  well  as  the  boundary  conditions  are  expressed  as  polynomials  in  e. 

The  discussion  of  the  parabolic  problem  [1]  gives  rise  to 
results  that  are  of  the  same  nature  as  in  the  elliptic  problem  [7] 
(except  that  in  this  case  Aronson  had  already  developed  a  method  for 
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showing  the  remainder  term  to  be  0(e),  which  in  this  thesis  is 

applied  to  the  elliptic  case) .  Once  again  a  method  is  devised  of 

improving  the  accuracy  of  the  representation  from  containing  an 

n 

error  term  which  is  0(e)  to  one  which  is  0(e)  as  e  -*  0  , 

and  this  type  of  analysis  extended  to  the  case  where  a,  b,  c,  f  and 

the  boundary  conditions  are  expressed  as  polynomials  in  e. 

The  conditions  under  which  the  results  of  the  elliptic 
problem  which  are  obtained  in  [7]  are  rather  restrictive,  and  so 
in  Chapter  IV  we  present  a  paper  by  Wasow  [10]  in  which  a  manner  of 
relaxing  these  restrictions  is  given  in  a  special  case.  In  this 
presentation,  we  append  a  comment  on  the  manner  in  which  the  sign 
of  e  as  e  0  determines  the  portion  of  the  boundary  over  which 
the  boundary  layer  phenomenon  will  be  exhibited.  We  also  show  how 
the  method  of  analysis  developed  by  Wasow  is  applicable  to  the  case 
of  the  corresponding  parabolic  problem. 


CHAPTER  II 


LEVINSON'S  METHOD  FOR  DEALING  WITH  THE  ELLIPTIC  SECOND  OPDER 
PDE  eAu  +  a(x,y)u  +  b(x,y)u  +  c(x,y)u  =  d(x,y) 

x  y 


1 .  Introduction 

In  1950,  Levinson  in  [7]  studied  the  asymptotic  expansion 
as  e  -*  0+  of  u(x,y;s),  the  solution  of  the  first  boundary  value 
problem  for 


eAu  +  a(x,y)ux  +  b(x,y)u^  +  c(x,y)u  =  d(x,y)  (1.1) 

for  small  e,  under  certain  mild  conditions  in  a  domain  D  .  The 
main  tools  of  the  proof  are  the  maximum  principle  for  elliptic 
differential  equations,  the  theory  of  characteristics  of  first  order 
partial  differential  equations,  and  an  integral  inequality  to  prove 
the  smallness  of  the  remainder.  His  error  or  remainder  estimate,  how¬ 
ever,  is  not  the  best  possible,  for  we  can  improve  it  from  0(/e’)  to 
0(e)  by  means  of  constructing  a  suitable  comparison  function  as  shown 
in  §7.  Furthermore,  in  §§8,9,  by  obtaining  more  terms  in  the  formal 

expansion  (5.2)  in  §5  and  (6.1)  in  §6  we  show  how  to  obtain  the 

N+l 

solution  with  an  error  term  equal  to  0(e  )  for  any  given  integer 

N  >  0  .  We  also  study  the  case  where  the  coefficients  of  the  elliptic 
equation  can  be  expanded  as  polynomials  in  e  . 

2.  Statement  of  the  Problem 


Let  D  be  an  open  simply  or  multiply  connected  region  in 


' 
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the  xy-plane.  F^  is  a  portion  of  the  boundary  of  D  which  is  nowhere 
parallel  to  the  characteristics  of  problem  below.  We  consider 

problems 


PDE  a(x,y)U  +  b(x,y)U  +  c(x,y)U  =  d(x,y)  (x,y)  e  D 
x  y 

IC  U  =  <f>  (x,y)  e  rQ 


(2.1) 


2* 


PDE  eAu  +  a(x,y)u  +  b(x,y)u  +  c(x,y)u  =  d(x,y) 


(x,y)  e  D  (2.2) 


BC 


<j>  (x,y)  £  r 

O 

$  (x,y)  e  QD-ro) 


where  ip  is  chosen  such  that  u  e  C(dD)  .  We  are  interested  in  the 
behavior  of  u(x,y;e)  as  e  -+  0+  and,  in  particular,  in  its  relation 
to  U(x,y)  of  . 

The  characteristics  curves ^  of  problem  P^  are  the 
solution  of 


dx  _  dy  , 

a(x,y)  “  b(x,y) 

First  we  require  the  following  hypotheses: 

H^.  Let  3D,  the  boundary  of  D,  consist  of  a  finite  number  of 
simply  closed  curves  3D^  (i  =  l,2,...,n).  Thus  IbF3D  is  a  closed 
connected  set.  Let  D+3D  be  contained  in  an  open  connected  region 


(1) 


See  [3] ,  pp .  62-64 . 


*S 
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6 

Dq  .  Assume  a(x,y),  b(x,y),  c(x,y)  and  d(x,y)  in  C  (Dq)  • 

3D^  (i  =  l,2,...,n)  are  given  by  functions  yielding  x  and  y  in 
terms  of  arc  length  along  each  closed  curve;  let  these  functions  and 
their  partial  derivatives  up  to  sixth  order  with  respect  to  arc  length 
be  continuous.  Let  the  boundary  values  assigned  to  u  along  each 
closed  curve  3D^  of  D  and  their  partial  derivatives  up  to  sixth 
order  with  respect  to  arc  length  be  continuous. 

2 

l^.  There  exists  a  function  A(x,y)  e  C  (Do)  such  that 


aA  +  bA  >0 

x  y 

2  2 

Remark  2.1.  implies  a  +  b 

2  2 

a  +  b  =0  at  some  point,  then  a 
too . 

^2  alternate»  c(x,y)  <  0  (x,y) 


(x,y)  e  Dq  . 

>  0  for  (x,y)  e  D  ,  since  if 

=  b  =  0  so  that  aA  +  bA  =  0 

x  y 

£  D 

o 


Remark  2.2.  Throughout  we  may  replace  by  alternate,  as 

theorem  3.2  shows. 


3 .  Maximum  Principle 

2 

Theorem  3.1.  Let  u(x,y)  e  C  (D)  satisfy 

Au  +  a(x,y)u  +  8(x,y)u  +  y(x,y)u  =  0 

x  y 

where  a,8,y  £  C(D)  ,  y  <  0  for  (x,y)  e  D  and  D  is  an  open 
domain.  Then  u  cannot  have  a  positive  maximum  nor  a  negative 


minimum  in  D  . 
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Corollary .  If  in  D 


Au  +  a(x,y)ux  +  8(x,y)uy  +  y(x,y)u  =  p(x,y) 

and  y_<-K<0,  |  p  |  <_  M  for  (x,y)  e  D  ,  then  u  cannot  have  a 

postive  maximum  nor  negative  minimum  in  D  exceeding  M/K  in 

magnitude.  If  u  e  C(D+8D)  and  max  |u(x,y) |  <_  M/K  ,  then 

(x,y)e8D 

max  | u(x,y)  |  <_  M/K  , 

(x,y)eD 


Proof.  If  u  has  a  positive  maximum  at  a  point  p  e  D  ,  then 

Au(p)  £  0  and  u  (p)  ®  u  (p)  »  0  .  Thus 

x  y 


(Au  +  au  4*  8u 
x 

Y (p)  •  u(p)  * 


+  yu)  *  p(p) 


p(p)  -  Au(p)  >_  p(p) 


Hence 


“(P)  1^7  i  M/K  • 


Similarly,  if  u  has  a  negative  minimum  at  the  point  p  ,  then 


Y (p)  *  u(p)  -  p(p)  -  Au(p) 
<  p(p)  • 


This  implies  that 


u(p) 


P(P) 

-y(p) 


>  -  M/k 


Given  and 


and 


u(x,y) 


e  C2(D) 


satisfies 


Theorem  3.2. 


■ 
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eAu  +  au  +  bu  +  cu  *  d 
x  y 

If  there  exists  a  constant  m  such  that 


max  |d(x,y)|  _<  m  and  max  |u(x,y)|  _<  m  , 

(x,y) e(D+3D)  (x,y)e3D 

then  there  exists  a  constant  K  independent  of  e  for  small  e  >  0 
such  that 


u(x,y)  |  <  Km  (x,y)  e  D 


A  A  (  x  y  ) 

Proof .  Let  v(x,y)  =  e  iy  *  u(x,y)  where  A  >  0  is  a  constant 
such  that 


By  h2 


A(aA  +  bA  )  >  0  . 

x  y 

there  exists  such  an  A  .  Then 


-AA 


u  =  e 


v  , 


so  that  v  satisfies 


eAv  +  (a-2eAA  )v  +  (b-2eAA  )v 
x  x  y  y 

+[c-A(aA  +  bA  )  +  eA^(A^  +  A^)  -  eAAA]v  =  d  *  e^ 
x  y  x  y 

Since  all  the  functions  in  the  coefficient  of  v  are  constant,  we  have 

that  for  a  sufficiently  small  e  the  coefficient  of  v  is  always 

negative  so  that  we  may  apply  the  corollary  in  this  section  to  get 

|v(x,y) |  is  bounded  in  D  .  Since  u  «  v  e  ^  ,  and  as  A  is  a 

-AA 

continuous  function  on  a  closed  bounded  domain,  e  is  bounded  above 
so  that  |u(x,y)|  is  also  bounded  above. 


4. 


Some  Conventions,  An  Existence  and  Uniqueness  Theorem 


Let  us  specify  the  positive  direction  on  3D  as  the  one 

taking  you  around  counter-clockwise  and  denote  the  arc  length  on 

3D  by  s  .  Let  T  be  an  arc  of  one  of  the  curves  of  3D  and  let 

o 

no  characteristic  of  (2.3)  be  tangent  to  .  Let  be  another 

arc  of  3D  such  that  the  characteristics  through  into  D 

intersect  and  so  that  also  nowhere  tangent  to  the  characteristic 

curves  of  (2.3).  Then  the  closed  region  bounded  by  T  ,  and  the 

two  characteristics  through  the  end  points  of  and  is  called 

a  "regular  quadrilateral". 


Since  neither  To  nor  has  any  characteristic  of  (2.3) 


tangent  to  it ,  we  have  on  3D 


b  i*  _  dy_  . 

b  ds  3  ds  *  U  ‘ 


If  on  T  we  have 
o 


bf- 

ds 


a^  <  0 

ds 


(4.1) 


then  on  r ^  we  have 


dx 

ds 


a  >  0 
ds 


(4.2) 


Theorem  4.1.  Given 

and 

and 

problem  of 

§2  .  Then 

there  exists  a  solution 

u(x,y;e) 

for 

any  small  e  >  0 

and 

u,u  ,u  e  C(IH-3D)  ,  u  e 
x  y 

c2(d)  . 
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Theorem  4.1  follows  from  results  of  Lichtenstein  [8]. 

Remark.  For  e  <  0  ,  the  roles  of  T  and  F,  are  interchanged. 

-  o  1 

For  a  proof  see  §2.2  ,  chapter  IV. 

5.  Statement  and  Proof  of  the  Main  Theorem  in  a  Special  Case 

Define  the  problem 

P_:  PDE  eAu  +  au  +  bu  +  cu  =  0  (x,y)  e  D  (5.1) 

j  x  y  j 

f0(s)  (x,y)  e  T- 

BC  u  -  < 

[  0  (x,y)  e  OD-r^ 

2 

whose  solution  is  u(x,y;s)  and  we  assumed  that  0(s)  e  C  (T^)  where 
s  is  arc  length  along  and  0(s)  =  0  over  a  small  segment  at  each 
end  of  ri  .  The  solution  of  the  associated  reduced  problem  is,  of 
course,  U(x,y)  =  0  . 

Main  Theorem.  Given  and  (or  alternate)  and  problem  P^. 

Then  in  any  regular  quadrilateral  Q  with  sides  and  ,  as 

above,  in  D+3D  we  have 

u(x,y ; e)  -  z(x,y;e)  +  eZ(x,y;e)  (5.2) 


where 


Z(x,y;e) 


f  0(1)  uniformly  in  Q  as  e  -+  0+ 
[_  0  on  3D 


and 
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-g(x>y) 

0 

z(x,y;e)  =  e  .  h(x,y)  on  and  at  points  in  Q  near  F^  (5.3) 


where 

(i)  g,h  e  C2(Q)  , 

(ii)  g  =  0,  h 

=  0(s) 

on  T±  ,  and 

(iii) 

g  >  0  inside  Q  . 

There  exists 

6  >  0 

(S  /  0 

such  that  z  =  0(e 

as  e 

->  0+  uniformly  in 

that  part  of  Q 

where 

(5.3)  no  longer  holds. 

Proof.  By  corollary,  §3  in  the  case  of  alternate  or  by  theorem  3.2 

in  the  case  of  ,  u(x,y;e)  is  uniformly  bounded  in  D+3D  as  e  ->  0"*"  . 

If  we  substitute  (5.2)  into  (5.1)  we  obtain 

7  Ig*  +  gy  -  Ugx  +  bgy)]  e"(g/e) 

+  [(a  -  2gx)hx  +  (b  -  2gy)hy  +  (c  -  Ag)h]e'(g/e) 

+  eAhe"^8^^  +  e[eAZ  +  aZ  +  bZ  +  c]  =  0  . 

x  y 

We  obtain  g(x,y)  by  setting  the  coefficient  of  e  equal  to 

zero  to  get  the  following  problem  for  g(x,y)  : 

2  2 

P,  :  PDE  g  +  g  -  (ag  +  bg  )  =  0  in  Q  near  ri 

4  ex  ey  &x  by  x  1 

IC  g  =  0  on  ri  and  g  >  0  in  Q  near  F^ 

—  (  p  /  0 

We  obtain  h(x,y)  by  setting  the  coefficient  of  e  equal  to 

zero  to  get  the  following  problem  for  h(x,y)  : 

P5  :  PDE  (a  -  2gx)hx  +  (b  -  2gx)hy  +  (c  -  Ag)h  -  0  in  Q  near  T± 

IC  h  -  0(s)  on  T1 

Now  we  solve  problem  P^  .  Let  P  =  gx  ,  q  =  g  .  Then  the 
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characteristic  system  for  is 


dx  _  - 

dt  ~  2p  ’ 

if  "  2q  '  b  > 

if  =  V  +  bxq  ’ 

if  ”  ayP +  y  > 

(5.4) 

-  P2  +  ,2  . 

2  2 

p  +  q  *  ap  +  bq 

We  shall  show  that  for  our  purpose  the  characteristics  of  (2.3) 
determine  the  segment  of  3D  over  which  we  make  use  of  (5.4)  and  that 
(4.2)  is  a  sufficient  condition  for  the  existence  of  the  solution  of 

P4  • 

To  get  g(x,y)  near  F^  in  Q  ,  it  is  sufficient  to  show 

that  no  characteristic  is  tangent  to  to  show  g  >  0  in  Q  near 

dfi  2  2 

ri  ,  observe  from  (5.4)  that  =  p  +  q  will  suffice  provided  p  =  0 

and  q  =  0  do  not  simultaneously  hold  and  dt  >  0  as  we  go  into  Q  . 

Now  on  T  we  have,  as  g  =  0  ,  that 

3g  .  _ 

if  =  pxs  +  qys  =  0  • 


From  (5.4)  we  have 

2  2 

p  +  q  =  ap  +  bq 
Thus  on  T±  we  may  solve  to  get 


x 
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Hence 


p  -  <ays  -  bxs)yg 

q  =  -(ay  -  bx  )x 
Js  s  s 


(5.5) 


(as  x  +  y  =  1  ,  and  we  reject  the  solution  p  =  q  =  0  ) .  By  (5.4), 
s  s 

(5.5)  and  (4.2)  we  obtain 


x  y_  -  xv  =  bx  -  ay  >  0 
st  f's  s  Js 


on 


Thus  no  characteristic  of  (5.4)  is  tangent  to  r  . 

implies  that  the  characteristics  of  (5.4)  go  into  Q 

2  2 

from  0  .  Again  (5.5)  implies  p  +  q  =  (ay  -  bx  ) 

s  s 

(Jcr 

Thus  >0  on  T1  and,  by  continuity,  in  Q  near 
g  =  0  on  T1  ,  g>0  in  Q  near  . 


Further 
as  t 
2  >  0 

ri  • 


(4.2) 
increases 
on  r1  . 
Hence,  as 


Let  the  arc  ,  ABCD  ,  have  A  and  D  as  its  end  points. 

ri  is  given  parametrically  as  x  =  x(s) ,  y  =  y(s)  .  0(s)  =0  on  AB 
and  CD  .  The  initial  values  at  t  =  0  on  for  (5.4)  read: 


x(s ,0)  =  x (s)  ,  y(s,0)  ■  y (s) 

g(s ,0)  =  0  , 

p(s ,0)  =  - [b (x(s) ,y (s)  •  x  (s)  -  ay  ]y  (s)  , 

S  SS 

q(s , 0)  =  [bx  -  ay  ]  *  x  . 

b  b  S 


Since  A(s,t)  =  ^  =  bx  -  ay  >0  on  r.  ,  then  by  continuity 

a  9Z )  S  S  X 

A(s,t)  >  0  for  some  region  extending  into  Q  .  As  a,  b,  x(s),  y(s) 
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are  all  of  class  by  ,  we  have  that  g  e  ^ 

Let  arcs  BF  and  CE  be  the  projections  of  the  characteristics 
of  (5.4)  into  the  xy-plane.  is  the  closed  region  ABCDEFA  ,  and 

AFED  is  sufficiently  close  to  that  lies  in  the  open  set  in 

which  g  exists.  Also  AFED  is  chosen  such  that 

(i)  x,y  along  AFED  expressed  in  terms  of  arc  length  are 

3 

of  class  C  , 

(ii)  the  normals  to  AFED  outward  to  are  directed 

into  Q  (this  is  arranged  by  making  [_  BAF  and 
l_  CDE  acute)  , 

2  2 

(iii)  AFED  is  sufficiently  close  to  that  p  +  q  >0 

2  2 

in  (this  is  possible  as  p  +  q  >0  on  ) 

(see  figure  2. 1) . 


(1) 


See  [2],  theorems  7.1  and  7.2,  pp.  24-25. 
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Since 

_ 

dt 

2  2 

p  +  q  >  0  , 

there  exists  6  > 

0  such  that 

g (x,y)  >  fit 

in  T 

and 

there  exists 

6  ^  >  0  such  that 

g(x,y)  _>  26x 

for  (x,y) 

on  arc 

EF  . 

is 


We  now  solve  problem  P 


5  ‘ 


The  characteristic  system  for 


dx 

dt 


-a  +  2gx  , 


dh 

dt 


(c-Ag)h 


iz 

dt 


-b  +  2g 


y 


(5.6) 


where  =  p  and  g^  ■  q  .  Since  the  first  two  equations  of  (5.6) 
and  the  first  two  equations  of  (5.4)  are  the  same,  the  characteristics 
of  P^  and  P^  in  the  xy-plane  are  the  same. 

The  last  equation  of  (5.6)  implies  that 

/*■  (c-Ag)dt 

h(s,t)  =  h(s)  e 

4  2  3 

Since  g  e  C  ,  then  Ag  e  C  .  Hence  h  e  C  for  t  ,  but 

2  2 

h  e  C  for  s  .  Therefore  h  e  C  .  h(s,0)  =  0(s)  =  0  for  s  on 

AB  and  CD  implies  that  h  =  0  in  the  curvilinear  triangles  ABF 

and  CDE  .  As  g,h  e  C  (T^)  then  z  =  h  *  e  &  e  C  (T^)  . 

z,  z  ,  z  ,  z  ,  z  and  z  vanish  on  AF  and  ED  and  are  all 
x  y  xx  xy  yy 

0(e  )  uniformly  on  EF  for  small  e  >  0  . 

We  now  extend  the  definition  of  z  a  little  further  into 
0  .  Let  a  be  the  arc  length  along  AFED  so  that  AFED  is  given 
parametrically  by  x  =  x(a)  ,  y  =  y(a)  .  We  erect  a  family  of  normals 
to  AFED  exterior  to  T^  with  x  as  arc  length  on  the  normals. 


' 
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These  normals  are  given  by 


x(t  ,o)  =  x(a)  -  y '  (or)  *  T  . 


y (t,o)  =  y (a)  +  x' (a)  •  x  . 


Since  ^ “  1  +  (x"y*  -  xfy")T  ,  on  x  =  0  (i.e.  AFED)  , 

d  \0  >  T  ) 

we  have  that  7- =1^0.  So  that  there  exists  a  1-1  correspondence 

9(a, t) 

between  (x,y)  and  (o,x)  near  t  =  0  .  Now  choose  such  that 

^  0  for  0  .1  T  .1  T-^  *  Let  GH  be  the  curve  along  which  we 

have  t  =  .  Let  be  the  curvilinear  quadrilateral  AFEDHGA  (Fig.  2.1), 

We  extend  z(x,y;e)  into  such  that  z  =  0  in  (Q  -  -  T^) 

2 

and  z  e  C  (Q)  as  follows. 


Define 


z  =  (1  ~  ^-)3  [z1(o)  +  t  •  z2(o)  +  t2  *  z3(o)] 


(5.7) 


2 

where  z..  (o)  =  z  ,  z  (o)  =  ~  +  —  z  (a)  and  z  (o)  =  -r-  -~j-  + 

1  Z  9  T  T  -  1  J  Z  ,,  Z 

1  dT 

3  3 

H - z_(o) - r  z.  (o)  on  t  =  0  .  The  functions  z,  z  ,  z  ,  z  ,  z  , 

t,2  21  x  y  xx  xy 

T 

1 

T 

z  all  vanish  on  the  arc  GH  ,  since  all  involve  a  factor  (1  -  — ) 

yy  t1 

Due  to  the  fact  that  z  =  0  in  the  curvilinear  triangles  ABF  and 


CDE  ,  z  and  its  partial  derivatives  up  to  second  order  vanish  on 

-  6l/e 


AF  and  DE  and  are  0(e 


)  for  (x,y)  e  (Q-T^  . 


Recall  that  by  (5.2)  we  have 


u(x,y ; e)  =  z(x,y;e)  +  eZ(z,y;e) 
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By  (5.3)  ,  z  =  e  6  •  h  in  where  g  and  h  are  solutions  of 

and  P  ,  respectively.  By  substituting  into  (5.1)  we  obtain 


eAZ  +  aZ  +  bZ  +  cZ 
x  y 


h  ‘ 


e(e/e) 


(x,y)  e  T1 


(5.8) 


Again  z  is  given  by  (5.7)  in  so  that  we  have 


eAZ  +  aZ  +  bZ  +  cZ  =  —  (eAz  +  az  +  bz  +  cz) 
x  y  e  x  y 


i  '  Vs 

=  0(±e  ) 


as  e  ■>  0+  (x,y)  e 


(5.9) 


Lastly,  z  =  0  in  (D+SD-T^-T^)  so  that  we  have 

eAZ  +  aZ  +  bZ  +  cZ  =  0  (x,y)  e  (Q-T.-T.)  .  (5.10) 

x  y  l  / 

Thus  from  (5.8),  (5.9)  and  (5.10)  we  have  that  Z  solves  the  problem 

PDE  eAZ  +  aZ^  +  bZ^  +  cZ  =  0(1)  uniformly  (x,y)  e  D+3D,  e  >  0 
BC  Z  =  0  (x,y)  e  3D 

Applying  the  maximum  principle  of  theorem  3.2  to  this  problem  we  have 
that 

Z (x,y ; e)  =  0(1)  (5.11) 

uniformly  in  D+3D  as  e  -*  0**"  .  Hence 

u(x,y ; e)  =  z(x,y;e)  +  0(e) 
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uniformly  in  D+3D  as  e  ->  0+  . 

6.  Statement  and  Proof  of  the  Main  Theorem  in  the  General  Case 

In  this  section,  instead  of  dealing  with  the  special  case 
of  problem  and  its  associated  reduced  problem  we  return  to  problems 
and  P^  • 

Main  Theorem.  Given  H^,  (or  alternate),  problems  P^  and 

P 2  with  solutions  U(x,y)  and  u(x,y;e),  respectively.  Then  in  any 
regular  quadrilateral  Q  with  sides  To  and  as  given  by  (4.1) 

and  (4.2)  in  D+3D  we  have  that 

u(x,y ; e)  =  U(x,y)  +  z(x,y;e)  +  eZ(x,y;e)  +  w(x,y;e)  (6.1) 

where  w  =  0(/e)  as  e  ->  0+  uniformly  in  Q  and  w  =  0  on 
and  r  . 

-g(x,y) 

z  =  e  £  •  h(x,y)  (6. 2) 

2 

on  T  and  at  points  in  Q  near  ,  where  g,h  e  C  (Q)  ;  g  =  0 

on  T  and  g  >  0  inside  Q  near  ;  h  =  ip  -  U  on  .  Z  =  0(1) 

uniformly  in  0  as  e  ->  0+  and  Z  =  0  on  3Q  .  There  exists  6  >  0 

_  6  f£ 

such  that  z  =  0(e  )  uniformly  in  that  part  of  Q  where  (6.2) 

no  longer  holds. 

Proof .  By  characteristics  in  this  section  we  mean  the  characteristics 
of  (2.3)  and  not  of  (5.4).  The  characteristics  of  (2.3)  can  be  designated 
by  the  coordinate  of  arc  length  s  on  .  Let  s  at  the  end  points 
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of  ri  be  given  by  s^  and  •  Since  (4.1)  and  (4.2)  are 
inequalities,  we  can  always  extend  r  and  thereby  enlarging 

Q  such  that  it  contains  the  characteristics  given 

s^  -  36  £  s  £  s2  +  36  for  some  6  >  0 


Denote  the  regular  quadrilateral  given  by  -  26  _<  s  _<  s^  +  26 
by  Q  and  the  segments  of  0  by  r  and  .  Let 

r  ^  -  U  on  r 
0(s)  =  /  1 

0  on  near  s^  -  36  and  s^  +  36 

2 

Let  0(s)  e  C  on  [s^-36,  s2+36]  .  Let  z(x,y;e)  and  Z(x,y;e) 
be  given  as  in  §5,  and  w(x,y;e)  be  given  by  (6.1).  w(x,y)  is 
uniformly  bounded  in  Q  as  e  0+  ,  since  u(x,y;e)  , 
z(x,y;e)  and  Z(x,y;e)  are  uniformly  bounded  in  Q  as 
Now  U(x,y)  satisfies 

(x,y)  e  Q 

E  Y 

O 

.  Also  z  is  given  by 

(x,y)  e  T1 

(x,y)  e  T2 

(x,y)  e  (Q-T1-T2) 


PDE  aU  +  bU  +  cU  =  d 

x  y 


IC  U  -  <J> 


(x,y) 


And  u(x,y;e)  satisfies  problem  P2 


e‘(g/e).h 


z  =  (  (6.7) 

0 


U(x,y)  , 


0 


where  g,h  solve  P^  and  P^  of  §5  ,  respectively.  Lastly  the  problem 
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Z  solves  is  given  in  §5  .  Hence  w  satisfies  the  problem: 


P,  :  PDE  eAw  +  aw  +  bw  4-  cw  =  -eAU  (x,y)  e  Q 


x 


(6.3) 


BC  w  =  0  (x,y)  e  (r  +  r^) 


w  in  general  is  non-zero  on  3D  -  (T  +  r.)  ;  indeed  w  =  u  there 

o  1 

2  — 

except  for  a  small  segment  at  each  end.  w  is  of  class  C  (Q)  ,  since 
z  e  C2(Q)  . 


Instead  of  using  the  characteristic  integral  inequality  of 
Levinson  [7]  to  obtain  that  w  =  0(/e)  in  Q  ,  we  shall  show  that 
in  next  section,  w  =  0(e)  in  Q  by  means  of  a  suitable  comparison 
function. 


7 .  The  Estimate  of  w(x,y) 
Let  Q 


0q  be  embeded  in  a  series  of  regular  quadrilaterals 


(k  *  0,1, 2, 3)  which  are  bounded  by  the  characteristics  s  ®  s^  - 


k6 


and  s  °  S2  +  k6  ,  and  the  sides 


(i) 


(i  =  1,2) .  Let  o 


(i) 


(i  =  1,2)  be  the  characteristics  bounding  (k  =  0,1,2)  (see 

figure  2.2). 
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Lemma  7.1.  Let  and  (or  alternate)  hold.  If 

W  =  W(x,y;e)  is  the  solution  of 

PDE  eAW  +  aW  +  bW  +  cW  =  0  in  0, 

x 

f W  =  0  on 
BC  / 

[  |  W  |  <  K  e11 

where  K  is  a  constant  independent  of  e  ,  then 

.  1 

n  +  2  + 

W  =  0(e  )  uninformly  in  as  e  ->  0  ,  (7.4) 

Proof .  Suppose  V  =  V(x,y;e)  solves  the  problem 


PDE  eAV  +  aV  + 

X 

bV  + 

y 

cV  =  0 

in  Qk 

(7.5) 

r  V  _>  0 

on 

r^1)  and 
k 

(7.6) 

BC  J 

\  V  >_  K  en 

on 

„  _  rt(l) 

o  -  ok 

and  o  »  a<2> 

k 

(7.7) 

1 

n  +  -R 

with  V  =  0(e 

)  uniformly  in 

j 

Qk  ^  as  e  ->  0 

(7.8) 

By  theorem  3.1,  V  >_  0 

in 

Qk  .  Now 

consider  the  function 

P(x,y;e) 

=  V  (x 

,y;e)  -  W( 

x,y;e)  . 

Then  P  solves 

PDE  eAP  +  aP  + 

bP 

+  cP  =  0 

in  Q. 

X 

y 

k 

BC  P  ^  0 

on  8Qk 

(1) 


and 


(2) 


on  a  =  a 


k 

(1) 


and 


(7.2) 


o  =  o. 


(2) 
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By  theorem  3.1,  P  >_  0  uniformly  in  .  Therefore  V  _>  W  uniformly 
in  .  Similarly,  by  considering  P  =  V  +  W  ,  we  find  that 
V  +  W  >_  0  or  W_>-V  uniformly  in  so  that  |w|  <_  | V j  .  Provided 

we  can  exhibit  a  V  such  that  (7.5),  (7.6),  (7.7)  and  (7.8)  hold, 
then 


n  +  i 

W  =  0(e  )  uniformly  in 


To  find  such  a  V  ,  consider  the  function 


V(x,y;e)  =  l 


2  -y(a) 

Tk  (a)  — 

e  •  n  (x,y;c)  +  /e  Y(x,y;s)  , 

a=l 


,  (a)  1  f  (a), 2 

where  Yk  =  —  (a-®k  ^ 

/e 


Substitute  into  (7.5)  to  get 


^  _y(a) 

\  e  k  [ ah^ J  +  bh^ J  +  { c+4 (o-a^7 )  (a^  +  c^)}hvu'7] 

a=l 


k  r_.(a)  ^  uu(a)  x  ^4.A^_a(a))(a2  +  a2xlt_(a) 

k  x  y‘ 


+  /eG(x,y;e)  4-  /e(£AY  +  aY  +  bY  +  cY)  =  0 

X  y  '  * 


(7.9) 


where 


G(x,y ; e)  = 


2 

l 

a=l 


■Y 


(a) 

k  r  t~  (a)  ,  ,  (a)w_  ,  (a)  .  _  ,  (a) 


e  [/s  An  -4(a-oy  )  (a  n  +  a  h  ) 

k  xx  y  y 


-2{ (a-a ) Aa  +  (a2  +  a2)}h^a^] 
k  x  y 


(7.10) 


(a)  COi') 

We  take  h  =  h  (x,y;e)  to  be  the  solution  of  the  initial  value 


problem 
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PDE 


,  (a)  ,  ,,  (a)  ,  r  .  .  ,  (a)w  2  .  2.  n,(a) 
ah  +  bh  +  [c+4  (c-a,  )  (a  +  a  )]h 

x  y  k  x  y 


0  in  Qk  (7.11) 


IC 


(a)  K  e 


n 


N 


on  r 


(1) 


where  N  >  0  is  a  constant  specified  in  (7.13).  Now  the  characteristic 
equations  of  (7.11)  are 


dx 

dt 


=  a 


dZ  = 

dt 


b  , 


dh 


(a) 


r  ,  /  /  (ot)  w  2  2.  ,  (a) 

^  =  [c  +  4(o-ok  )(ax  +  oy)]h 


(7.12) 


Hence  the  characteristics  of  (7.11)  on  the  xy-plane  are  the  same  as 

(a) 

those  of  the  reduced  equation  (2.1).  Hence  h  exists  uniquely 
in  and  is  in  class  C^(Qk)  .  By  (7.12),  and  t  =  0  on  , 

,  ,  „  /  t[c+4(a-a^a))(o2+o2)]dt 

(a)  _  K  n  o  k  xy 

h  -  N  £ 


f  \  /  \  T*  f  cdt 

_  (a)  ,  (a)  _  K  n  o 

On  o  -  o,  ,h  =  —  e  e 
k  N 


.  Thus  if  we  let 


N  =  min  e 
(a »t)  e  Q, 


/  cdt 
o 


(7.13) 


then  0  <  N  <  1  and  h^  >  K  en  on  a  =  o^a^ 
—  —  k 


n. 


From  (7.10),  G  e  C  (Q^)  and  G  =  0(e  )  uniformly  in  . 


From  (7.9),  Y  solves 


PDE  eAY  +  aY  +  bY  +  cY  =  -G  in  Q. 

x  y  k 


BC 


Y  =  0 


on  9Q 


k  ’ 
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Now  by  theorem  3.1,  Y  =  0(en)  uniformly  in  for  small  enough 

e  >  0  .  Thus  we  have  found  the  required  V  . 

Theorem  7.1.  Let  H-^  and  (or  alternate)  hold,  and 

w(x,y;e)  be  the  solution  of 

PDE  eAw  +  aw  +  bw  +  cw  =  -eAU  (x,y)  e 

x  y  y  2 

/  w  =  0  on  and 

w|  <  B  on  a  =  and  a  = 

where  B  is  independent  of  e  .  Then  w  =  0(e)  uniformly  in  Qq 

n+ 

as  e  ->  0 

Proof .  Since  we  are  interested  only  in  the  order  of  magnitude  of 
w  ,  we  can  assume  without  loss  of  generality  that  c(x,y)  <  0  in 
Q2  .  For  if  it  is  not,  we  need  merely  use  the  substitution  of 
theorem  3.2  to  insure  it. 

Now  let  Wq  =  Wo(x,y;e)  solve  the  problem 

PDE  eAW  +  aW  +  bW  +  cW  =  -eAU  in  Q. 

o  o,x  o,y  o  x2 

BD  W  =  0  on  3Q„ 

o  2 

By  theorem  3.2,  Wq  =  0(e)  uniformly  in  .  Put  w  =  Wq  +  W  . 

Then  W  solves  the  problem 
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PDE  eAW  +  aW  +  bW  +  cW  =  0 

x  y 


in  Q. 


BC 


on  and 


on  a 


(1) 


and  o=o 


(2) 


Then  by  lemma  7.1  with  n  =  0  and  k  =  2  we  have  that  W  =  0(/e) 
uniformly  in  0^  ,  i.e.  there  exists  K  >  0  ,  independent  of  e 
such  that  | W J  <  K  /e  in  .  By  applying  lemma  7.1  again  with 
n  =  1/ 2  and  k  =  1  ,  it  follows  that  W  =  0(e)  uniformly  in  0^ 
as  e  ->  0+  . 


8.  A  Refinement  of  the  Special  Case  of  §5 


In  section  5  we  have  a  solution  of  problem  consisting 

of  a  boundary  layer  term  plus  an  error  or  remainder  term  eZ  which 
is  0(e)  in  D+3D  as  e  ■*  0+  .  In  this  section  we  shall  improve  the 
order  of  the  remainder  from  1  to  N+l  in  e  as  e  -*  0+  .  The  notation 
of  §5  will  be  retained,  except  for  slight  modifications  introduced  at 
the  appropriate  place.  In  particular,  hypothesis  is  replaced 

by  hypothesis  Hj  in  which  we  have: 


(i)  a,b,c,d  e  C2N+6  (Dq) 

(ii)  The  functions  yielding  x  and  y  in  terms  of  arc 

length  specifying  each  closed  curve  are  now  of  class 

2N+6  ,  -i  .. 

C  in  arc  length. 

QXT  I  o 

(iii)  6(s)  e  C  (T^)  anc^  t^e  boundarv  condition  for  u 
is  of  class  C2^+2  (3D)  . 
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Theorem  8.1.  Given  H|,  (or  alternate)  and  problem  with 

6(s)  modified  as  indicated  in  hypothesis  .  Then  in  any  regular 

quaduilateral  Q,  with  sides  T  and  T,  as  before,  in  D+3D  we 

o  1 

have  that 


u(x,y;e)  -  z(x,y;e)  +  eN+1  Z(x,y;e) 

where  Z  =  0(1)  uniformly  in  Q  as  e  -*  0+  and  Z 
3D  ,  and 


vanishes  on 


(8.1) 


z 


_  g.(x,£)  N 

£  I  en  hn(x,y) 


(8.2) 


on  r±  and  at  points  in  Q  near  ,  where 

(a)  g  e  C5+2N(Q)  ,  hn  e  C3+2(N_n)(Q)  (n  =  1,2,...,N) 

(b)  g  =  0  ,  h  =0,  h  =0  (n  =  1,...,N)  on  T.  , 

on  i 

(c)  g  >  0  inside  Q  . 

J  iS  >  0,  such  that  z  =  0(e  )  as  e  -*  0  uniformly  in  that 

part  of  Q  where  (8.2)  no  longer  holds. 


Proof.  In  the  region  T^  we  have  z  now  given  by  (8.2).  Thus  the 
substitution  of  (8.1)  into  (5.1)  will  yield 


±  e'(g/e)[(g*+gy)-(agx+bg  )]  l  £nhn 

n=o 

,  ,  v  N  3h  3h 

+  e'  (g  £  l  en[(a-2gx)  +  (b-2g  )  +  (c-Ag)hJ 

n=o  y 

+  e”^8^  l  en+1  Ah  +  eN+1[eAZ+aZ  +bZ  +cZ]  =  0  . 

L  n  x  y 

n=o 
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As  before,  we  obtain  g(x,y)  as  the  solution  to  problem  and  hQ(x,y), 
which  is  the  same  as  h(x,y),  as  the  solution  to  problem  in  §5  . 

For  hn(x,y)  we  examine  the  coefficient  of  en  to  obtain  problem 


P 


* 

n 


3h 

PDE  (a-2g  )  — tt-  +  (b-2g  )  +  (c-Ag)h  =  -Ah  1  in  Q  near  T. 

&x  8x  &y  n  n-1  1 

IC  h  =  0  on  T.  . 

n  1 


Now  as 


a,b,  x (s ) ,  y(s) 


e  C 


6+2N 


we  have  that 


g  e 


c5+2N 


(1) 


The  characteristic  LaGrange  systems  for  problems  P  and  P^ 
are  respectively 


/  dx  - 
dt  “  2gx 


dt  =  2gy  '  b 


s  dx 

r  aT  ■  2gx  -  a 


■  2gy  '  b 


dh 

_£  .  (c-ag)ho 


dh 


n 


dt  =  (c-A8>hn  +  AVl 


Hence  the  characteristics  in  the  xy-plane  are  the  same  for  all  these 
problems.  Since 

dh 

-3!  ■  (c_Ag)ho 


,  r5+2N 

and  g  e  C 

,  then 

x  _3+2N  ..  .  ,  4+2N 

Ag  e  C  ,  so  that  h  e  C 

o 

in  t  and 

,  3+2N 

h  e  C  in 

o 

s  .  A 

similar  argument  applies  to  h^  . 

Thus 

h  Ec3+2(N-n) 

n 

• 

^See  [2]  pp. 

24-25. 
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Now 


/  (c-Ag)dt 

hQ(t,s)  *  hQ(0,s)  *  e 


But  h  (0,s)  *  6(s)  =  0  for  s  on  AB  and  CD  .  Thus  h  =  0  on 
o  o 

the  curvilinear  triangles  ABF  and  CDE  .  Again  the  problem 


dh 


T~-  05  (c-Ag)h  +  Ah  i 

dt  °  n  n-1 


ODE 


h 


IC 


=  0 


n 


has  solution 


/Qt  (c-Ag)dt  rt  /Qt  (c-Ag)dt 


hn(t,s)  -  e 


Ah  -  (t,s)dt 
n-1 


e 


But  as  Ah  _  =  0  in  the  curvilinear  triangles  ABF  and  CDE  , 


so  is  h  =  0  there, 
n 

The  extension  of  z  into  T^  in  §5  remains  essentially 
the  same.  However  the  estimation  of  Z  is  changed  to  the  following. 
Recall  that  z  is  defined  as  follows: 


(x,y)  e  T1 


0 


(x,y)  e  (D+aD-T^T^ 


Thus  the  PDE  for  Z  is  given  by 


■ 
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r  -Ahn  •  e  (g/°  (x,y)  e  T1 

I  l  e-(g/£)  + 

eAZ+aZ  +bZ+cZ  »  ( — r-  (eAz+az  +bz  +cz)  =  0( - as  e  -*  0  in  T0 

x  y  \  N+l  x  y  N+l  2 

e  J  e 

^  0  (x,y)  e  (D+3D-T;l-T2) 

Provided  e  >  0  is  sufficiently  small,  the  R.H.S.  is  0(1)  in  D+3D  . 
Hence  we  can  write  the  problem  for  Z  as  follows : 

PDE  eAZ  +  aZ  +  bZ  +  cZ  =  0(1)  uniformly  in  D+3D,  e  >  0 

x  y 

BC  Z  =  0  (x,y)  e  3D  . 

By  the  maximum  principle  of  theorem  3.2,  we  have  that 

Z(x,y;e)  =  0(1)  as  e  -*  0+ 
uniformly  in  D+3D  . 


9 .  A  Refinement  of  the  General  Case  of  §6 

In  section  6,  we  have  a  solution  of  problem  p in  a  regular 
quadrilateral  Q  of  the  form: 

u(x,y ; e)  -  U(x,y)  +  z(x,y;e)  +  eZ(x,y;e)  +  w(x,y;e)  (9.1) 

where  (i)  U  is  the  solution  of  the  reduced  equation, 

(ii)  z  is  the  boundary  layer  term  and  in  the  region  is 

j  u  -(g/e). 

given  by  z  »  e  h  , 

(iii)  Z  is  the  remainder  term  associated  with  this  form  of  the 
boundary  layer  and  is  0(1)  uniformly  in  D+3D  , 


■ 
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(iv)  w  is  the  remainder  term  introduced  by  the  non-homogeneous 

term  d  in  the  equation  (2.1)  and  the  initial  condition 

U  =  <p  on  T  ,  and  is  0(e)  . 
o 


In  this  section  we  increase  the  order  of  our  approximation  from  1 
to  N+l  in  e  as  e  ■+  0+  in  a  regular  quadrilateral  0  by  assuming 
a  solution  of  form 


N 

u(x,y ;  e)  =  UQ(x,y)  +  l  e11  Un(x,y)+z(x,y;e)+eN+1Z(x,y;e)+eNw  (9.2) 


n=l 


where  (i)  u  is  the  solution  to  problem  under  the  hypothesis 
rather  than  , 

(ii)  z,  Z,  w  play  the  same  role  as  in  §6  although  they  have 
different  formulations  here, 


(iii)  U  is  the  same  as  U  in  §6  , 
o 

N  n 

(iv)  Z,  e  U  is  a  term  whose  introduction  is  necessary  in 
n=i  n 

order  that  the  above  terms  have  the  properties  stated. 
(More  precisely  the  counter  part  of  (6.3)  above,  which 
is  (9.5)  below,  requires  a  term  of  order  N+l  in  e  on 
the  R.H.S.) 


Theorem  9.1.  Given  (or  1^  alternate)  and  problems  and 

P ^  with  solutions  UQ(x,y)  and  u(x,y;e)  ,  respectively.  Then  in  any 
regular  quadrilateral  0  with  sides  and  as  given  by  (4.1) 

and  (4.2)  in  D+3D  we  have  that  w  is  given  by  (9.2)  where 
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(i)  Z  =  0(1)  uniformly  in  D+ 3D  and  Z  =  0  on  3D  . 

(ii)  w  =  0(e)  uniformly  in  Q  as  e  ->  0+  and  w  =  0  on 

T  and  r,  . 
o  1 

_  N 

(iii)  z  =  e  E  l  A  (x,y)  on  ri  and  at  points 

n=o 

in  Q  near  (9.3) 

where  (a)  g  e  C  (Q)  ,  h  e  C  (0), 

n 

(b)  g  =  0  ,  h  -  ip-u,  h  =  -U  (n  =  1 ,2 , . . .  ,N) 

o  n  n 

on  T1  , 


(c)  g  >  0  inside  Q  near  T. 


(iv)  U  e  C2+2(N-n)(Q)  . 
n 


-  «1/c  + 

6^  >  0  such  that  z  =  0(e  )  as  e  ■>  0  uniformly  in  that 


part  of  Q  where  (9.3)  no  longer  holds. 


Proof .  We  shall  use  the  notation  of  §6  again  except  where  noted. 
The  substitution  of  (9.2)  into  (2.2)  yields 


N  N  3U  3U 

e  [  enAU  +  l  en(a  ~  +  b  -r^  +  c  U  ) 
L  n  L  3x  3x  n 

n=o  n=o 


N 


+  “  e  ^^[(g^  +  gy)  -  (agv  +  bgTr)]  l  enh 


x  y  °  n 

n=o 


+  e‘(8/E>  l  en  [(a-2gx)  ^  +  (b-2g  )  +  (c-Ag)hnl 

n=o 


+  e'(g/e)  l  en+1  h  +  eN+1[eAZ+aZ  +bZ  +cZ] 
L  n  x  y 

n=o  J 


‘ 


, 
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N 

+  e  [eAw  +  aw  4-  bw  +  cw]  =  d 
x  y 


(9.4) 


As  before  Uo  solves  the  reduced  problem  P^  (and  is  thus  the  same 
as  U  in  §6).  Un(x,y)  solves  problem  P”  (n  =  1,2,...,N)  : 


3U 


PDE 


n 

3x 


3U 


+  b 


+  cU 

3y  n 


-AU 


n-1 


in  Q+BQ 


IC 


U  -  0 
n 


on  r 


o 


The  function  g(x,y)  is  again  the  same  as  in  §5  and  hQ(x,y)  agrees 


» ?  t 


with  h(x,y)  of  §5  .  And  the  functions  h  (x,y)  solve  problem  P 

n  n 

(n  =  1,2, ..  .  ,11)  : 


ah 


ah 


PDE  (a-2g  )  — +  (b-2g  )  — +  (c-Ag)h  =  -Ah  ..  in  0  near  T 
x  3x  y  ay  n  n-1  1 


IC 


h  =  -U 
n  n 


on  r. 


for  the  function  Z(x,y;e)  we  have  the  problem 


-AhN  e 


-(g/e) 


in  T, 


PDE  eAZ+aZ  +bZ  +cZ 
x  y 


[eAz+azx+bzy+cz]  in  T  2 
in  (Q+aq)  -  (Tj+T^ 


BC 


Z  =  0 


on  aq 


Note  that  while  it  is  natural  to  set  Z  =  0  on  T  and  we 

o  1 

arbitrarily  set  Z  =  0  on  the  balance  of  c>0  .  Further,  by  choosing 
g  as  given  in  §5  sufficiently  small,  we  can  insure  that  the  R.H.S. 

jL 

of  the  PDE  is  0(1)  .  So  that  by  the  maximum  principle  Z  =  0(1)  in  0  . 
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Lastly,  as  u,  E  e*1  U  ,  z,  Z  are  all  defined,  so  is  w 

n=o  n 

as  the  solution  of  the  problem 

PDE  eAw  +  aw  +  bw  +  cw  =  -eAU..  in  Q+30  (9.5) 

x  y  N 

BC  w  —  0  on  (T  +  T  )  (9.5) ' 

o  1 

w  in  general  is  non-zero  on  30  -  (T  +  T,)  .  Since  this  is  the  same 

o  1 

as  (6.3)  and  (6.3)’  with  here  replacing  U  ,  the  argument  for 

w  in  §7  may  be  repeated  to  show  that  w  =  0(e)  in  (0+3Q)  . 


10.  An  Extension  of  the  Special  Case  of  §5 

We  wish  to  extend  our  problem  to  the  case  where 


M 


i(x,y;e)  =  J  e  am(x,y) 


m=o 


and 

b,  c, 

d  are  similarly  defined,  as 

are  the  boundary  conditions 

for 

u  on 

3D  .  The  problem  P^  for 

u(x,y ;  e) 

is 

PDE 

eAu  +  au  +  bu  +  cu  =  0 
x  y 

(x,y) 

e  D  (10.1) 

r  9 (s ; e) 

(x,y)  e 

k 

V 

rl 

BC 

^  0 

(x,y)  e 

OD  -  T±) 

where  we  assume  that 


N 


(i)  9  (s ;  e)  =  l 


n=o 


en  0  (s)  with  0  e  n^(r.) 

n  n  1 


' 
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(ii)  0  -  0  over  a  small  segment  at  each  end  of 

Note  that  the  associated  reduced  problem  has  the  solution 
We  set 

u(x,y;e)  =  z(x,y;e)  +  eN+1  Z(x,y;e) 

N 

where  z  =  e  ^(x>y)/£)  J  e11  h  (x,y)  with  N  >_  M 

n=o 


r 


l  ' 


U(x,y;e)  =  0  . 


(10.2) 


M 


Z-  I 


em  Z  (x,y;e) 
m 


m=o 


Substitute  (10.2)  into  (10.1)  to  get 


(g/e)  j  en+1[^h  -  —  Ag  *h  -  — (g  h  +g  h  )  +  At  (g2+g2)h  ] 

L  n  e  &  n  e  x  n,x  by  n,y  2  °x  y  n 

n=o  e 


,  ,  .  W-N  m 

+  e  j  I  em(a  h  +b  h  +c  h  ) 

“  “  m-n  n,x  m-n  n,y  m-n  n 


m=o  n=o 


M+N  m 


_  e  £  J  em  1(a  g  +b  g  )h 

L  L  m-n  x  m-n  y  n 


m=o  n=o 


+  e 


N+2  "r  m 


M  M . i  2M  m 

J  emAZ  +eN'  J  j  em(a  Z  +b  Z  +c  Z  ) 

L  m  L  L  m_n  n  x  m_n  n  y  m_n  n 

m=o  m=o  n=o 


=  0 


g(x,y)  solves 


PDE 


g2+g2-(ag  +  b  g  )  =  0 

°x  6y  o6x  o  y 


IC 


g  =  0  on  r. 
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g  >  0  in  0  near  .  Hence  g  is  the  same  as  in  problem  . 

The  coefficient  of  en  e  ^ ^  ,  for  n  =  0,1,..., N  yields  the 


following  equation  for  h  : 

n 


(a  -2g  )h  +  (b  -2g  )h  +  (c  -ang  -b-g  -Ag)h 
o  °x  n,x  o  &x  n,y  o  l&x  l6y  6  n 


n 


_Ah  .  -  J  (a.  h  .  +  b,  h  .  +  c,  h  ) 

n-1  k  n-k,x  k  n-k,y  k  n-k 


n 


+  l  (ak+l8x  +  bk+l8y)hn-k  ’ 
k=l  J 


with  initial  conditions  h  =  6  on  T.  .  The  coefficient  of  e 

n  n  1 


N+l 


is 


eAZ+aZ  +  b  Z  +  c  Z 
o  oo,x  oo,y  oo 


M-l 


M 


£  ^ak+l8x+bk+l8y')hN+l-k  J  (akhN+l-k,x+bkhN+l-k,y 
k=l  k=l 


+CkhN+l-k')  AhN 


r(§/£) 


(10.3) 


The  coefficient  of  eN+^+m  for  m  =  1  M-l,  yields  for  Z  : 

m 


eAZ+aZ  +  b  Z  +  c  Z  =  e  (g/e) 
m  o  m,x  o  m,y  o  m 


M-l 

^ak+l8x+bk+l8y^bN+l+m-k 


M 

1^=1  N+l+m-k,x  bkbN+l+m-k,y  +  ScbN+l+m-k,y^ 


n 


^  ^ak^m-k,x  +  bk^m-k,y  +  Ck^m-k^ 


k=l 


(10.4) 
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Lastly,  the  partial  differential  equation  for  is 


M 


eAZv  +  a  +  b  Z„  +  c  Z„ 

M  o  M,x  o  M,y  o  M 


k^![akZM-k,x  +  bkZM-k,y  +  ckZM_k] 


M  M 

y  em  y  [a,  Z  .  +  b.  Z  .  +  c.  Zx,^_  ,  ]  . 

,L  k  M+m-k,x  k  M+m-k,y  k  M-hn-k,y 
m=l  k-m 


(10.5) 


Now  in  the  region  we  define 


N 


:<o,T;e)  =  l  £n  Zn(a,T!e) 


n=o 


where 


zn(a,T;e) 


(1  - 


3  2 

77>  I 


m=o 


m  f  v 

x  z  (a;e) 
m  n 


— (g/ £ ) 

continues  the  solution  e  h  from  T-  into  T_  such  that 

n  1  2 

2 

z  e  C  (T,  +  T-)  .  Thus  the  remainder  function  Z  satisfies  the 
12  o 

PDE  (10.3)  in  the  region  ,  (10.6)  below  in  the  region  T ^  ,  and 

(10.7)  in  the  region  D  -  (T^  +  T^)  where  we  have 


eAZ  +  a  Z  +  b  Z 
o  o  o,x  o  o,y 


—  —  i  [eAz  +  az  +  bz  +  czl 
N+l  x  y 


=  0 


-(6^/ e) 
N+I 


as  e  ■>  0 


(10.6) 


and 


eAZ  +  a  Z  +  b  Z  +cZ  =  0 
o  oo,x  oo,y  oo 


(10.7) 
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with  the  boundary  condition  Zq  =  0  on  3D  .  Hence  we  have  that 
Zq  =  0(1)  as  e  ->  0+  uniformly  in  the  region  D  . 

In  similar  fashion,  we  have  that  Z  (m  =  1,2,...,M-1) 

m 

satisifes  (10.4)  in  the  region  while  in  the  region  D  -  we 

have 

m 

eAZ  +  a  Z  +  b  Z  +  c  Z  -  -  J  (a.  Z  ,  +b.  Z  .  +c.  Z  .  ) 

m  o  m,x  o  m,y  o  m  ,  k  m-k,x  k  m-k,y  k  m-k 


Since  the  boundary  conditions  are  Z^  =  0  on  3D  ,  we  have,  by 
induction,  that  Z^  =  0(1)  as  e  -*  0+  uniformly  in  the  region  D  . 


Finally,  Z^  satisfies  (10.5)  in  the  region  D  so  that, 

with  the  boundary  condition  Zw  =0  on  3D  ,  we  have  the  Zw  =  0(1) 

M  M 

as  e  0  uniformly  in  the  region  D  .  Thus  we  conclude  that  the 


remainder  term  Z  is  0(1)  uniformly  in  D  as  e  -*  0  .  As  before, 

,  ,  v  t  v  6+2N  .  5+2N 

we  assume  a,  b,  c,  x(s),  y(s)  e  C  so  that  g  e  C  , 

,  3+2 (N-n)  .  n  ,  M  ,  _  r2  ,  .  3+2 (N-n) 

h  e  C  for  n  =  0,1,..., N  ,  and  Z  e  C  where  0  e  C 

n  n 


11 .  An  Extension  of  the  General  Case  of  §6 

As  in  §10,  we  wish  to  have  the  coefficients  a,  b,  c,  d 
expressed  as  polynomials  of  order  M  in  terms  of  the  parameter  e  . 

We  shall  then  construct  a  solution  of  form 

N 

u(x,y;e)  =  U  (x,y)  +  l  e  un(x,y)+z (x,y ; e)  +  e  Z(x,y;e)  +  e  w(x,y;e 

n=l 

The  discussion  is  almost  the  same  as  in  §9  except  now  U  solves 

n 


(n  =  1,2,... ,N) 
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n 

PDE  a  U  +  b  U  +  c  U  =  -AU  ,  -  J  (a,  U  , 

o  n,x  o  n,y  o  n  n-1  k  n-k,x 

+  bkVk,y  +  ckUn-k>  ln  «*«» 

IC  U  =  0  on  r 

n  o- 

For  h  ,  the  initial  condition  becomes  h  =  9  -  U  on  r,  . 

n  n  n  n  1 

Throughout  the  entire  region  D  we  add  to  the  R.H.S.  of  equations 
(10.3)  and  (10.4)  that  is  the  coefficient  of  e^4™  (m  =  1,2,..., M) 
the  term^ 

M 

l  (a!  U\t i  i  +  b.  U.^  .  +  c.  UXT .  .  )  . 

.  L  k  N+m-k,x  k  N+m-k,y  k  N+m-k 
k=m 

Lastly,  as  the  problem  solved  by  w  in  §9  is  unchanged,  we  have  once 
again  the  w  =  0(e)  uniformly  in  D  as  e  ->  0+  . 


(1) 


Note  that  this  does  not  change  the  order  of  (m  =  0,1, . . . ,M-1) . 


CHAPTER  III 


LINEAR  PARABOLIC  DIFFERENTIAL  EQUATION 
CONTAINING  A  SMALL  PARAMETER 


1 .  Introduction 

Aronson  in  [1]  investigates  the  first  boundary  value  problem 
for  the  parabolic  equation 

eu  +  a(x,y)u  -  b(x,y)u  +  c(x,y)u  =  d(x,y)  (1.1) 

XX  x  y 

in  D,  an  open  region  in  xy-plane  bounded  above  by  the  line  y  =  y2 

and  below  by  the  concave-upward  curve  T,  where  e  >  0  is  a 

parameter  and  b(x,y)  >  0.  If  the  given  data  is  sufficiently 

smooth,  there  exists  a  unique  solution  u(x,y;e)  satisfying  equation 

(1.1)  and  the  boundary  condition  u  =  4>  on  T.  It  is  shown  that 

the  solution  u(x,y;e)  has  the  form  u(x,y;e)  =  U(x,y)  +  z(x,y;e)  +  w(x,y;e) 

in  a  regular  multilateral  M  (which  will  be  defined  later)  where  U 

is  the  solution  of 

a(x,y)Ux  -  b(x,y)U^  +  c(x,y)U  «  d(x,y)  (x,y)  e  D  (1.2) 
u  =  c()  (x,y)  e  T  . 

w  is  0(e)  in  M  and  z  is  a  boundary  layer  term  similar  to  that 
obtained  by  Levinson  [7]. 


(1) 


See  [4] ,  p .  41 . 


' 
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In  §§7,  8,  9  and  10,  we  shall  extend  Aronson’s  results 
as  we  have  extended,  in  the  last  chapter,  Levinson’s  results. 

2 .  Statement  of  the  Problem 

Let  Dq  be  an  open  region  in  the  xy-plane,  and  let 
and  r3  be  two  closed  arcs  in  which  are  defined  for  <_  y  <_ 

and  do  not  meet  for  y  >  y^.  Denote  by  DC  D^"^  the  closed  region 
bounded  above  by  a  segment  1^  of  y  =  y^,  below  by  a  segment 
of  y  -  y^,  and  laterally  by  T ^  and  (See  Fig. 3*1).  Let  s 

denote  arc  length  on  3D  =  4-  1^  measured  positively 

in  the  counter-clockwise  direction.  We  consider  the  problem 


P  • 
1  * 


PDE 

eu  +  au  - 

XX  X 

bu 

y 

+  cu  =  d 

IC 

( 

-  (j>^ 

on 

h 

u  =  , 

*2 

on 

~2 

-  *3 

on 

T3 

<*V 

4,  £  c4 

in  s 

and 

*9  "  V 

(x,y)  £  (D  +  3D) 


r2nrx,  r3n  ri  ,  respectively, 


(1) 


If  A  is  an  open  set,  then  by  A  we  mean  the  closure  of  A. 
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Let  us  also  assume  the  following  and  call  them  hypothesis  H^: 


(A) 

a,  b,  c 

4 

,  d  e  C 

<v 

(B) 

e  >  0, 

b(x,y)  > 

m  >  0  in  Dq 

(C) 

T  2  and 

~3  are 

4 

of  class  C  with  respect  to  s  and 

are  nowhere  horizontal. 


The  characteristics  of  the  reduced  partial  differential 
equation  (1.2)  are  given  by 


dx 

dt 


a(x,y) 


= 

dt 


b(x,y)  . 


(2.1) 


In  view  of  parts  (A)  and  (B) ,  there  exists  a  unique  solution 

of  (2.1)  through  every  point  of  which  can  be  continued  as 

far  as  the  boundary  of  and  which  is  continuous  with  respect  to 

its  initial  values. 


Let  I  be  a  closed  sub-arc  of  r  (i  =  1,2,3)  such  that 
no  characteristic  is  tangent  to  I  and  such  that  every  characteristic 
which  starts  on  I  enters  D  for  increasing  t,  and  first  leaves 

D  via  a  closed  subarc  J  C  +  T_  +  T.  -I.  Let  no  characteristic 

2  3  4 

starting  on  I  be  tangent  to  J.  We  call  the  closed  region  bounded 
by  I,  J  and  the  characteristics  joining  their  end-points  a  regular 
multi-lateral  M  (see  Fig.  3.2). 


If  BD  is  given  parametrically  as  x  =  x(s),  y  =  y(s),  then 


A 


(2.2) 


. 
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is  (i)  equal  to  0  at  any  point  T  is  tangent  to  a  characteristic 
of  (1.2),  (ii)  is  less  than  zero  on  I,  and  (iii)  is  greater  than 
zero  on  J.  Thus  every  point  in  D  +  r  which  does  not  lie  on  a 
characteristic  which  is  tangent  to  T  or  on  the  characteristics 
through  O  T  (i  =  2,3)  can  be  included  in  a  regular  multilateral. 

Main  Theorem.  Let  hold.  Consider  the  problem 

Pn  PDE  aU  -  bU  +  cU  =  d  (x,y)  e  M  (2.3) 

2  x  y 

IC  U  =  u  (x,y)  £  I 

and  problem  P^  above.  Then  in  the  multilateral  M, 

u(x,y ; e)  =  U(x,y)  +  z2(x,y;e)  +  z3(x,y;e)  +  w(x,y;e)  (2.4) 

where  (i)  w  =  0(e)  uniformly  in  M  as  e  ->•  0+,  w  =  0  on 
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m  n  {v±  +  r2  +  r3>, 

(ii)  if  J  n  T  =  0  (i  =  2,3),  then  z  =  0;  if 
J  fl  T  f  0  (i  =  2,3),  then 

g±Cx,y) 

zi(x,y;e)  =  e  £  •  h^x.y)  (2.5) 

at  point  of  M  on  or  near  JOT.  ,  where  g.  =  0  on  J  f|  r., 

i  i  l 

-  -  2 
g,>0  off  J  n  r . ,  h.=u-U  on  J  n  r .  and  g.,  h.  e  C  . 
l  i  i  i  i  i 

C  6  /  s  ^ 

Moreover,  there  exists  6  >  0  such  that  z^  =  0  (e  )  uniformly 

in  that  part  of  M  where  (2.5)  no  longer  holds,  and  z^  =  0  on  I. 

3 .  Maximum  Principles 

Let  D  be  the  region  given  at  the  beginning  of  §2,  and 

consider 

L  [u]  =  eu  +  a(x,y)u  +  3(x,y)u  +  y(x,y)u  =  6(x,y)  (3.1) 

G  XX  X  y 

in  D.  u(x,y;e)  is  said  to  be  a  regular  solution  of  L£(u)  =  6 

in  D  +  8D  if  u  satisfies  the  differential  equation  in  D  and  is 

continuous  in  D  +  3D  and  if  u  ,  u  ,  u  e  C(D  +  T.). 

x  y  xx  4 

Lemma  3.1.  Assume 

(I)  a,  3,  y,  <5  e  C1 (II) (Ill)(Dq)  . 

(II)  e  >  0  and  3  >  m  >  0  in  D^. 

(Ill)  and  are  piecewise  smooth  and  nowhere  horizontal. 

Functions  d> .  e  are  given  on  T.  (i  =  1,2,3)  such 

l  l 

that  <p  =  <j>^  at  r  D  (i  =  2,3). 


(IV) 
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Then,  for  fixed  e,  (3.1)  has  a  regular  solution  in  D  +  3D, 
which  is  equal  to  <j)  on  F  (i  =  1,2,3). 

Lemma  3.2.  Let  u  be  a  regular  solution  of  Le[u]  =0  in 
D  +  3D,  where  £  >  0,  3  >  0,  y<0  in  D+3D,  and  a,  3,  y 
are  bounded  in  D  +  3D.  Then  u  cannot  take  on  a  positive  maximum 
or  a  negative  minimum  in  D  +  r^. 

The  proofs  of  Lemmas  3.1  and  3.2  are  given  in  [4],  Chapter  2. 


Lemma  3.3.  Let  u(x,y;e)  be  a  regular  solution  of  (3.1),  where 
£  >  0,  3  >  0  in  D+3D,  and  a,  3,  Y  are  bounded  in  D  +  3D.  If 

1 6  (x,y)  |  <_n  in  D  +  3D  and  |  u  |  <_  n  on  r^+^+T^,  the*1 
there  exists  a  constant  K  >  0,  independent  of  e,  such  that 
|u|<_Kn  in  D  +  3D. 

lev 

Proof .  Let  u  =  e  y  v(x,y;e)  where  k  >  0  is  a  constant  such 
that  y(x,y)  -  k  •  3(x,y)  <_-y<0  in  D+3D  for  some  constant 
y  >  0.  Then  (3.1)  becomes 


£ v  +  av  -  Sv  +  (y  -  k3)v  =  6  •  e 
xx  x  y 


-ky 


(3.2) 


Now  v  may  not  have  a  positive  maximum  nor  negative  minimum  in 

n  "kyl 

D  +  such  that  |v|  >  —  e  .  For  suppose  it  did  have  a  positive 

^  -ky- 

n  1 

maximum  at  P  £  D  such  that  v  >  —  e  .  Then 
o  y 


v  (P  )  <  0,  v  (P  )  =  v  (P  )  =  0. 
xx  o  —  x  o  y  o 


And 
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[ev  +  (y  -  kB)v] 

XX 


n  'kyl  'kyl 

<  -  \i  •  —  •  e  =  -  ne 


o 


-ky 


1 


i.e.  the  L.H.S.  of  (3.2)  is  less  than  -  ne  .  But  the  R.H.S. 


of  (3.2)  is  greater  than  or  equal  to  -  ne 


“ky. 


-ky, 


Similarly  for  a  negative  minimum,  v  <  -  —  e 


n  ^  1 


y 


on 

h  + 

r  +  r 

2  3* 

1  v  1 

<  n 

-kyi 

.  e 

k(y-yv 

1  u  1 

<  n 

•  e 

4. 

The 

Statement 

Now  |  u |  <_  n  on  implies  |v|  <_  ne 

-kyx 

in  D  +  F. ,  then 
4 


-ky, 


_n 

y 


l- 


max  {1,  — }  in  D  +  T,  +  r.  +  r_  +  I\.  Therefore 
,  y  14  2  3 

•  max  (1,  — }  <  k  •  n  in  (D  +  BD) . 

y  - 


Let  J *  CZ  r  (i  =  2,3)  be  any  closed  sub-arc  such  that 


A  >  0  on  J* 


(4.1) 


where  A  is  given  by  (2.2).  Let  0(s)  e  C  (J*)  such  that  0(s)  =  0 
in  a  neighborhood  of  each  end-point  of  J*.  Consider  the  problem 


P3: 

PDE 

eu 

XX 

+  au 

X 

bu 

y 

+  cu 

IC 

u  = 

r  0  (s) 

on 

J* 

l  o 

on 

h  ' 

The 

associated 

reduced 

equation 

of 

(4.2) 

(4.2) 


condition  has  solution  U(x,y)  =0  in  D  4-  3D. 

Theorem  4.1.  Given  and  problem  P^.  Then  in  any  regular 


multilateral  M  there  exists  a  regular  solution  of  form 


u(x,y;e)  =  z(x,y;e)  +  e  Z(x,y;e) 


(4.3) 


where 


(i)  Z(x,y;e) 


4  0(1)  uniformly  in  M  as  e  -+  0+ 
\  o  on  r1  +  r2  +  r3 


_  $(.x^.y) 

(ii)  z  =  e  £  h(x,y) 


(4.4) 


on  J* 

and 

at 

points  in  M 

near 

J*  where 

g,  h  e  C2(M) ; 

•k 

o 

n 

t»0 

h  = 

e 

on  J*  and 

g  >  0 

inside 

M 

near  J*.  There 

exists 

6  > 

0 

such  that  z 

=  0  (e 

-(6/e)} 

as 

e  ->  0+  uniformly  in 

that  part  of  M  where  (4.4)  no  longer  holds. 


Proof .  The  proof  paraphases  the  proof  in  §5,  Chapter  II  except 
for  the  appropriate  minor  changes  noted  below. 

By  lemma  3.3,  u(x,y;e)  is  uniformly  bounded  in  M  +  3M 
as  e  -*■  0+.  If  we  substitute  (4.3)  into  (4.2)  we  obtain 

7  [§x2  "  a8x  +  b§y]  e  ^g/E)  +  t(a  -  2§x>hx  "  bhY  +  (c  -  gxx)h]  e  (g/ 

+  eh  +  e[e  Z  +  aZ  -  bZ  +  cZ]  =  0. 

xx  xx  x  y 

We  obtain  g(x,y)  by  setting  the  coefficient  of  e  equal 

to  zero  to  get 

2 

p,  .  PDE  g  -  ag  +  bg  =0  in  M  near  J* 

4  •  °x  x  °y 

IC  g  =  0  on  J*  and  g  >  0  in  M  near  J*. 
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We  obtain  h(x,y)  by  setting  the  coefficient  of  e  equal  to 

zero  to  get 

Pr  •  PDE  (a  -  2g  )h  -  bh  +  (c  -  g  )h  =  0  in  M  near  J* 

5  °x  x  y  °xx 

IC  h  =  0  on  J*. 


The  method  of  solution  for  g(x,y)  of  and  h(x,y)  of  P^ 

is  the  same  manner  as  that  in  §5,  Chapter  II.  But  instead  of  (5.4), 
Chapter  II,  we  have  that  the  system  of  characteristic  equations  for 


P,  is 
4 


dx  n 

dt  =  2p  ~  a’ 


&  = 


dt 


b  , 


&  =  n2 

dt 


P  » 


dp 

■jr  =  a  p  -  b  q  , 
dt  xr  x 


and  we  obtain 


d  q  ,  2 

jT  =  a  p  -  b  q,  P  =  ap  -  bq  ; 

u  l  y  y 


p  =  (ays  +  bxs)/ys 


q  =  -x  (ay  +  bx  )/y 
n  s  s  s  s 


Instead  of  (5.6),  Chapter  II,  we  have  here  that  the  characteristic 
system  for  P^  is 


dx  ,  n 

— —  =  -a  +  2q  , 
dt  Hx  ’ 


dt 


b  , 


dh  ,  Nu 

—  =  (C  -  gxx)h. 


The  equations  that  Z  satisfies  in  various  regions  in  D  +  9D  are  the 
same  except  for  the  change  of  the  operator  from  elliptic  to  parabolic. 
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Also  throughout  the  whole  proof  we  replace  Q  and  F^  by  M  and  J*  , 
respectively. 


5.  The  Proof  of  Main  Theorem  in  the  General  Case 


Let  M  be  a  fixed  regular  multilateral  and  let  J  be  given 
in  terms  of  arc  length  for  s  1.  s2*  Since  each  characteristic 

of  (1.2)  in  M  meets  J  exactly  once,  the  coordinate  s  denotes  a 
unique  characteristic  in  M.  Since  (2.2)  involves  strict  inequalities, 
and  since  no  characteristic  in  M  meets  3D  -  (I  +  J),  we  can  always 
extend  I  and  J  to  obtain  an  enlarged  regular  multilateral  M^ 
which  contains  the  characteristics  s^  -  36  <_  s  <_  S2  +  36  for  some 
6  >  0.  Thus  M  =  Mq  can  always  be  embedded  in  a  series  of  regular 
multilaterals  (k  =  0,1, 2, 3)  which  contain  the  characteristics 

s,  -  k6  <  s  <  s~  +  k6.  Let  I.  and  J,  denote  the  non-characteristic 
boundaries  of  M  .  Let  U(x,y)  be  the  solution  of  which  takes  on 


the 


values  of  u(x,y;e)  on  I^.  Then  U  exists  uniquely  in  and 


is  of  class  C  (M2^  * 


Now  we  define  w  in  (2.4)  and  distinguish  the  following 


three  cases: 


Case  1.  J  c  T. .  Then  we  define  w(x,y;e)  by  u(x,y;e)  =  U(x,y)  + 
-  4 

w(x,y;e)  where  u  is  the  unique  regular  solution  of  problem 
P  in  D  +  3D. 

Case  2.  J  c  T.  (i  =  2,3).  Let  9(s)  =  u  -  U  on  J 2  and  let  0  be 

2 

continuous  over  -  J2  such  that  0  e  C  (J^)  and  0=0 
in  a  neighborhood  of  the  end-points  of  J^.  Let  u(x,y;e)  be 
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Case  3. 


the  regular  solution  as  obtained  in  §4  to  the  problem 


PDE  eu  +  au  -  bu  +  cu  =  0  in  D  +  3D 
xx  x  y 


IC 


0  on  J. 


u  = 


0  on  +  r2  +  r3  ”  * 


We  define  w  in  by 


u(x,y;e)  =  U(x,y)  +  u(x,y;e)  +  w(x,y;e). 

_  _ 

J  fl  T,  ^  0,  J  (f  r..  Let  D  c  D  be  the  closed  region  bounded 
4  4  o 

below  by  above  by  a  segment  of  the  line  y  =  y^  where 

— * 

^3  >  ^2*  anc*  laterally  by  the  arcs  r  ( i  =  2,3)  which 
are  obtained  by  continuing  the  F  above  y  =  y9  in  such  a 
way  that  hypothesis  H^(C)  continues  to  hold.  Now  H  /  0 

(k  =  0,1, 2, 3)  for  i=2  or  i  =  3  or  i  =  2  and  3.  For 

_  ^  _ 

such  an  i  we  have  J„  H  r .  c  V .  and  A  >  0  on  H  r.. 

2  i  i  2  i 

*  ; k 

By  continuity,  there  exists  a  closed  subarc  c;  r_^  such  that 

J.  coincides  with  J _  fl  r.  below  T.  fl  r.  ,  which  extends 

l  3  i  4  i 

—  —  * 

a  finite  distance  above  F.  fl  F.  and  A  >  0  on  J.  .  Let 

4  i  l 

0^(s)  =  u  -  U  on  fl  ri,  and  let  0^  be  continued  over 

*  —  2  *v 

J ,  -  such  that  0^  e  C  (J\)  and  0^  =  0  in  a  neighborhood 

* 

of  each  end-point  of  Let  u_^  =  u_^(x,y;e)  be  the  regular 

solution  of  the  problem 


PDE 


IC 


eu  4-  au 
xx  x 


* 

bu^  +  cu  =  0  (x,y)  e  (D  +  3D  ) 


u  = 


■0 . 

i 


0 


* 


(x,y)  e  J_ 


— *  — *  * 

(x,y)  s  (r1  +  r2  +  r3  -  j.). 
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^  ^  j*  jj 

Note  that,  as  (D  +  3D  )  c:  D0  and  A  >  0  on  JL  c  r  , 

the  results  of  §4  can  be  applied  to  u..  If  J.  H  r.  =  0, 

i  k  i 

we  set  u_^  =  0.  We  now  define  w  in  by 


u(x,y;e)  =  U(x,y)  +  u2(x,y;e)  +  u3(x,y;e)  +  w(x,y;e) . 


These  cases  are  shown  in  the  Figure  3.3  below.  The  arc  abc  is  a 
characteristic  of  (1.2)  which  is  tangent  to  at  b.  Note  that  in 

n 


Figure  3,3 

each  case  M  can  be  embedded  in  a  family  of  regular  multilaterals 
of  the  same  type. 


In  all  the  above  cases,  w  solves  the  problem 


PDE  ew  +  aw  -  bw  +  cw  =  -  eU  (x,y)  e  M_ 

xx  x  y  xx  2 

ic  w  =  o  (x,y)  e  m2  n  [r^  +  r2  +  r3] . 

w  is  unspecified  but  bounded  on  the  boundaries  of  M2  that  are 
characteristics,  w  is  uniformly  bounded  in  M?  as  e  ->  0,  since 
u,  u2,  u3  are  all  uniformly  bounded  in  D  +  3D  as  e  ->  0  and 


U  is  also  bounded. 


. 
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In  the  next  section  we  shall  show  that  w  =  0(e)  uniformly 
in  =  M  as  e  ->  0+.  This,  together  with  Theorem  4.1  in  §4  will 
complete  the  proof  of  the  main  theorem  in  §2. 


6 .  The  Proof  that  w  =  0(e) 

The  proof  in  this  section  is  analogous  to  §7,  Chapter  II. 
The  partial  differential  equation 


a(x,y)  o  -  b(x,y)  a  =  0 
x  y 


(6.1) 


has  the  solution  a(x,y)  e  C  (t^)  and  a  =  constant  are  the 
characteristics  of  (1.2).  Let  a  =  and  a  =  be  the 


characteristics  bounding  (k  =  0,1,2).  Let  y  be  the  largest 

* 

constant  such  that  the  line  y  =  y  has  at  least  one  point  in 

A 

common  with  ^ .  Let  be  the  closed  region  obtained  by 

*  — 

excluding  the  part  of  above  the  line  y  =  y  .  If  fl  f  0, 

A  _  A 

then  =  M^.  If  fl  =  0,  then  c  c  M^.  These  two 

cases  are  shown  in  Figure  3.4  and  Figure  3.5  respectively. 


Figure  3.4 


Figure  3.5 


X 
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Lemma  6.1.  Let  hold.  If  W(x,y;e)  is  a  regular  solution  of 


PDE 


IC 


eW  +  aW  -  bW  +  cW  =  0  in  M, 
xx  x  y  k 


f  w  =  o  on  n  [r1  +  r2  +  r3] 


|  I W I  <  Ken  on  a  =  and  o  =  o in  M 

^  k  k  k 

A 

where  K  is  independent  of  e  and  c(x,y)  <  0  in  M^,  then 

W  =  0(en+^)  uniformly  in  ^  as  £  -*0+, 

~  ~  * 

Proof .  Let  denote  the  boundary  of  below  the  line  y  =  y  . 

Suppose  V  =  V(x,y;e)  is  a  regular  solution  of  the  problem 


PDE 


IC 


with 


eV 

+  aV 

-  bV 

+ 

cV  =  0  in 

XX 

X 

y 

V  1 

0 

A 

on 

n  [rx  +  72 

V  > 

Ke11 

on 

o  = 

(D 

a,  and. 

k 

(2) 

’k  in  \ 


(6.2) 


V  =  0(en+lj)  uniformly  in  as  e  0+. 


By  lemma  3.2,  V  >_  0  in  M^.  Now  consider  the  function  P(x,y;e)  =  V  -  W. 
Then  P  solves 


PDE 


eP  +  aP  -  bP  +  cP  =  0  in  M. 
xx  x  y  tc 


IC  P  >_  0  on  3M^. 

A 

By  lemma  3.2,  P  >_  0  uniformly  in  M^.  Therefore  V  >_  W  uniformly 

A 

in  M  .  Similarly  by  considering  P  =  V  +  W,  we  find  that  W  >_  -V 

K 

A 

uniformly  in  M^.  Thus  the  lemma  will  be  proved  once  a  suitable  V  is 
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exhibited . 


Consider  the  function 

2  -y(o) 

V(x,y;e)  -  l  e  k  h^“hx,y;e)  +  y'z  Y(x,y;e)  , 
a=l 


where  y 


(a)  _  1  /-  _(aK2 

k  -  7=  (0  -  °k  >  • 


2  -r 

l  e 

a=l 


(a) 

k 


[ah 


(a) 


x 


Substitute  V  into  (6.2)  to  get 


, ,  a  ,  r  .  ,,  (a) . 2  2,  (a), 

bh  +  (c  +  4  (a  -  a.  )  a  }  h  ] 
v  k  x 


+  /T  G(x,y;e)  +  (eY  +  aY  -  bY  +  cY)  =  0  , 

xx  x  y 


2  -Y(a) 

where  G(x,y;e)  =  T  e  ^  [/e"  n  -  4 (a  -  a1t'a',)  a  h' 

xx  k  x  x 

a=l 


(a)s  _  ,  (a) 


-  2{o2  +  (a  -  a,(a))a  )  h(a)]  . 

x  k  xx 


(ot ) 

We  take  h  to  be  the  solution  of 


PDE 


IC 


ah^a^  -  bh^0^  +  [c  +  4 (a 
x  y 


(a) v  2  2,  ,  (a) 

°k  )  °x]  h 


=  0  in 


(a)  Ke 


n 


N 


on  I, 


where  N  >  0  is  a  constant  to  be  specified  later.  As  in  §7,  Chapter 
IT,  0  <  N  _<  1,  h^>_  Ken  on  a  =  in  M^,  if  we  let 

ft  cdt 

J  o  XT 

N  =  min  „  e  .  Also  Y  =  0(e  )  uniformly  in  M,  .  Thus 

(a , t)  e 

we  have  found  the  required  V. 


Let  hold.  If  w(x,y;e) 


Theorem  6.1. 


is  a  regular  solution  of 


' 
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PDE 


ew  +  aw  -  bw  +  cw  =  -  eU  (x,y)  £  M~ 

xx  x  y  xx  2 


IC 


w  =  o  on  m2  n  [r  +  r2  +  r  ] 

|  w |  <  B  on  a  =  and  a  =  in 


where  B  is  independent  of  e,  then  w  =  0(e)  uniformly  in 
as  e  ->  0+. 


Proof ♦  By  lemma  3.1  there  exists  a  regular  solution  Wo(x,y;e)  of 

PDE  eW  +  aW  -  bW  +  cW  =  -  eU  in  M0 

o,xx  o,x  o,y  o  xx  2 

IC  W  =  0  on  9M0  . 

o  2 

By  lemma  3.3,  =  0(e)  uniformly  in  •  Let  w  =  +  W.  Then 

W  solves 

A 

PDE  eW  +  aW  -  bW  +  cW  =  0  in  M0 
xx  x  y  2 

(  w  =  o  on  m2  n  [r1  +  ?2  +  T3 ] 

IC  i  i  ,  (1)  (2)  A 

y  |  W  |  <  B  on  0  =  0^  and  a  =  a ^  in  . 

By  lemma  6.1  with  n  =  0  and  k  =  2  we  have  that  W  =  O(e^)  uniformly 
in  M^;  i.e.  there  exists  K  >  0,  independent  of  e,  such  that 
|  W  |  <  K  /e"  in  M  .  Since  c  ,  by  lemma  6.1  with  n  =  y  and 
k  =  1,  it  follows  that  W  =  0(e)  uniformly  in  as  e  -*  0+. 


7 .  A  Refinement  of  the  Special  Case  of  §4 

Just  as  in  the  case  of  the  elliptic  problem,  we  wish  to 
improve  the  order  of  the  remainder  term  eZ  from  1  to  N  +  1  in  e 
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as  e  ■>  0  .  To  do  so  we  must  refine  our  hypothesis  to  H|  : 

(i)  a,  b,  c,  d  e  C^+2N(D) . 

(ii)  The  functions,  yielding  x  and  y  in  terms  of  arc 

length  specifying  each  boundary  curve  are  now  of  class 

o4+2N  i  -u 

C  in  arc  length. 

2N  * 

(iii)  0(s)  e  C  (J  )  and  the  initial  condition  for  u  are 
also  of  class  C2^(T). 


Theorem  7.1.  Given  H|  and  problem  P^.  Then  in  any  regular 
multilateral  M,  there  exist  a  regular  solution  of  form 


u(x,y;e) 


/  \  N+l  ,  N 

z(x,y;e)  +  e  Z(x,y;e) 


(7.1) 


where 


Z(x,y;e) 


0(1)  uniformly  in  M  as  e  -*■  0+ 
o  on  r1  +  r2  +  r3 


and 

2  =  e-g(x,y)/e  f  c"hn(x>y)  (7.2) 

n=o 

*  *  3+2N 

on  J  and  at  points  of  M  near  J  where  g  e  C  (M)  , 

h  e  c1+2(N"n) (M) ;  g  =  0,  h  =0,  h  =0  (n  =  1 ,2 , . . . ,N)  on  J*; 
n  on 

* 

g  >  0  inside  M  near  J  .  Moreover, there  exists  6  >  0  such  that 
(S  / 

z  =  0(e  )  as  e  -*  0  uniformly  in  that  part  of  M  where  (7.2) 

no  longer  holds. 


Proof.  In  the  region  T^  we  have  that  z  is  given  by  (7.2)  .  Thus 
the  substitution  of  (7.1)  into  (4.2)  will  yield 
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i  »-s/£  ro2  - 


N 


[gx  ■  (agx  +  b8y)]  l  e”  hn 

n=o 


N 


+  e  g/E  l  en  [(a  -  2g  )  h  -  bh  +  (c  -  g  )  h  ] 

°x  n  X  \r  °  w  -n 


n=o 


n,y 


xx  n 


N 


+  e"g/e  y  en+1  h  +  £N+1[eZ  +  aZ  -  bZ  +  cZ]  =  0 

L  n,xx  xx  x  y 

n=o 


As  before  we  obtain  g(x,y)  as  the  solution  to  problem  and 

h  (x,y)  as  the  solution  to  problem  Pr.  For  h  (x,y)  we  examine  the 
o  5  n 

coefficient  of  e11  to  obtain  (n  =  1,2,...,N) 


P*  :  PDE  (a  -  2g  )h  -  bh  +  (c  -  g  )h  =  -h 

n  x  n,x  n,y  &xx  n  n-l,xx 

5 

in  M  near  J 


IC  h  =  0  on  J 
n 

Now  as  a,  b,  x(s)  ,  y(s)  z  q^+2N^  we  have  that  g  e  C^+^N.  Again, 

the  characteristics  in  the  xy-plane  for  P.  and  P'  (n  =  0,1,..., N) 

4  n 

are  the  same.  Also,  h^  z  £l+2(N  n\  The  partial  differential  equation 
solved  by  Z  is 


eZ  +  aZ  -  bZ  +  cZ 
xx  x  y 


h.7  e 
N,xx 


- (g/e) 


in  T. 


— ,T  [zz  +  az  -  bz  +  cz] 
,N+1  xx  x  y 


=  0( 


~(&/z) 

N+l 


)  as  e  ->  0  in  T, 


0  in  (D  +  3D)  -  (T  4-  T2) . 


Thus  the  problem  for  Z  may  be  expressed  as 
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PDE  eZ  +  aZ  -  bZ  +  cZ  =  0(1)  uniformly  in  (D  +  3D) 
xx  x  y 

as  e  -*■  0+ 


ic  z  =  o  on  r1  +  r2  +  r3  . 


By  the  maximum  principle  of  lemma  3.3,  Z(x,y;e)  =  0(1)  as  e  -*  0+ 
uniformly  in  D  +  3D. 


8 .  A  Refinement  of  the  General  Case  of  §5 

The  modification  in  this  case  is  essentially  the  same  as 
given  in  §9  for  the  elliptic  case.  We  shall  have  that 


N 


u(x,y;e)  =  UQ(x,y)  +  l  en  U^(x,y)  +  u2(x,y;e)  +  u3(x,y;e)  +  eN  w(x,y;e) 


n=l 


where  u2  and  u3  are  defined  as  in  §5  but  with  a  remainder  term 
e^+^Z  obtained  as  in  §7. 


As  before,  Uq  solves  the  reduced  problem  P2  (and  is  thus 
the  same  as  U  in  §5)  while  Un(x,y)  solves  problem  P^  (n  =  1,2,...,N) 


PDE 

aU 

bU 

+  cU  =  -  U  , 

in 

(M  + 

3M) 

n,x 

n,y 

n  n-l,xx 

IC 

U  =  0 
n 

on 

I. 

(  QJ  ) 

The  function  g  J  ( 

x,y) 

associated  with  u 

a 

(a  = 

1,2) 

is 

again  the 

same  as 

in  §4,  and 

h(c) 

0 

(x,y)  agrees  with 

h(a) 

(x,y) 

of 

§4.  The 

function 

h^U.y) 

solves  the  problem  P”' 

(n  = 

1,2,. 

•  •  j 

N) 

PDE  (a  -  2g  )h(a)  -  bh(a)  +  (c  -  g  )h(a) 
°x  n,x  n,y  xx  n 


(a)  * 

-  hv  t  in  M  near  J,  v 
n-l,xx  (a) 
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IC 


(a) 


n 


-  U  on  J ,  s 
n  (a) 


(oi) 

For  the  function  Z  we  have  the  problem 


_hM  e-(g  M 

N  ,xx 


in  T, 


PDE  eZ(a)  +  aZ(a)  -  bZ(a)  +  cZ^ 
xx  x  v 


4T[^(0‘)+a^“)-bz<“)+cii 

N+l  xx  x  y 


V  0  in  (M  +  3M)  -  (T  +  T  ) 


in  T, 


IC  Z  =  0  on  ri  +  r2  +  r3  ‘ 


By  choosing  6^  as  given  in  §4  sufficiently  small,  we  can  insure  that 
the  R.H.S.  of  the  PDE  is  0(1)  so  that,  by  the  maximum  principle 
of  lemma  3.3,  Z  =  0(1)  in  M  . 

N 

r»  XI 

As  u,  )  e  U  ,  u^,  u_  are  all  defined,  so  also  is  w 

L  n  2  3 

n=o 

as  the  solution  of 


PDE  ew  +  aw  -  bw  +  cw  =  -  elL.  in  (M  +  3M) 

xx  x  y  N,xx 


IC 


w  =  0  on  I  +  J 


w  in  general  is  non-zero  on  3M  -  (I  +  J  ) . 


Since  this  is  the  same  as  that  of  §6,  we  once  again  obtain  that  w  =  0(e) 
uniformly  in^  M  as  e  ->  0  . 


9 .  An  Extension  of  the  Special  Case  of  §4 


We  wish  to  extend  our  problem  to  the  case  where 
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M 


a(x,y;e)  =  £ 


em  am(x,y) 


m=o 


and  b,  c,  d  are  similarly  defined,  as  are  the  boundary  conditions 
for  u  on  3D. 

The  problem  for  u(x,y;e)  is 


PDE 


eu  +  a  u  -  b  u  +  cu  =  0 
xx  x  y 


(x,y)  e  D 


IC 


u  = 


6 (s;e) 


0 


(x,y)  e  J  ) 
(x,y)  e  (3D 


* 

J  ) 


(9.1) 


where  we  assume  that 


N 


(i)  0 (s;e)  -  l  en  0  (s)  with  0  e  C1+2(N  n) ( J*) 

u  n  n 

n=o 


(ii)  0=0  over  a  small  segment  at  each  end  of  J  . 

Its  associated  reduced  problem  has  the  solution  U(x,y;e)  =  0  in 
D  +  3D. 

We  shall  construct  u(x,y;e)  of  the  form 


u(x,y;e)  =  z(x,y;e)  +  eN+1  Z(x,y;e) 


(9.2) 


where 


N 


_  e  (g(x,y)/e)  ^  £n  j^(x>y)  with  N  >_  M 


n=o 


M 


Z  "  l 


m=o 


m  _  .  v 

E  zm  (x.y;0 


and 


60  - 


Substituting  (9.2)  into  (9.1)  yields  the  following  results  in  the  region 


T 

L1  * 


g(x,y)  solves 


PDE  g  -ag  +  b  g  =0 
x  ox  o  y 

IC  g  =  0  on  J  and  g  >  0  in  M  near  J 


Thus  g  is  the  same  as  in  problem  P^.  h^  solves 


PDE  (a  -  2g  )h  -  b  h  +  (c  -  g  -  a.  g  +  b-g  )h 
o  °x  n,x  o  n,y  o  6xx  l&x  1  y  n 

n 


h  -  T  (a,  h  .  -  b.  h  .  +  c.  h  ) 

n-l,xx  k  n-k,x  k  n-k,y  k  n-k 


n 


+  E  ^ak+l8x  "  bk+l8y)  hn-k  * 
k=l  J 


IC  h  ■  6  on  J 
n  n 


And  the  partial  differential  equations: 

eZ  +  a  Z  -bZ  +  c  Z  =  e"('g/e') 
o,xx  oo, x  oo, y  oo 

M 


M-l 

^  (ak+l8x  "  bk+lSv')hN+l-k 


k-1 


"  ^^A+l-k^  +  Wl-k^  +  °khN+l-k)  “  \ , 


xx 


(9.3) 


eZ  +  a  Z  -  b  Z  +cZ*e 
m,xx  om,x  om,y  om 


“(g/e) 


M-l 


l  <al,4.l8„  -  ^1,4.1  8„)h 


k=l 


k+l°x  k+l6y  N+l-hn-k 


M 


I  ^akbN+l+m-k,x  hkhN+l+m-k , y  +  CkbN+H-m-k) 
k=l 


M 


-  J  (a,  Z  .  -b,Z  .  +  c.  Z  ,) 

L  k  m-k,x  k  m-k,y  k  m-k 

lv=  X 


(9.4) 
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M 


£ZM,xx  +  aoZM,x  boZM,y  +  CqZM  " 


-  I 


k-1 


^  akZM-k,x 


b.  Z  .  +  c  Z  } 

k  M-k,y  k  M-k 


M  M 

V  m 

L  e 

m=l  k=m 


J  {a,  Z  .  -  b.  Z _  .  +  C,  Z  ,  } 

L  k  M+m-k , x  k  M+m-k , y  k  M+m-k 


(9.5) 


For  the  region  ,  we  define 


N 


(o,t;e)  =  l  enz  (a,x;e) 
“  n 


n-o 


x  3 

where  z  (a,x;e)  =  (1 - )  T  x  z  (a;e)  continues  the  solution 

n  x  '  L  mn 

1  m-o 

-(e/e)  2 

e  6  h  from  T1  into  T.  such  that  z  e  C  (T.  +  T-). 

n  l  2  l  z 


The  remainder  Z  satisfies  the  PDE  (9.3)  in  T,  (9.6) 

o  i 

below  in  T^,  and  (9.7)  in  D  -  (T^  +  T^)  where 


eZ  +  a  Z  -bZ  +CZ 
o,xx  o  o,x  o  o,y  o  o 


-~b  {ez  +•  az  -  bz  +  cz} 
,N+1  xx  x  y 


0 


-(61/e) 

i 

N+l 


as  e  ■>  0 


(9.6) 


eZ  +  a  Z  -bZ  +  c  Z  =  0 
o,xx  o  o,x  o  o,y  oo 


(9.7) 


With  the  initial  condition  Z  *  0  on  T,  we  have  that  Z  =  0(1) 

o  o 

as  e  -*  0+  uniformly  in  D.  In  the  similar  fashion,  we  have  that 

Z  (m  =  1,2,...,M-1)  satisfies  (9.4)  in  T.  while  in  D  -  T.  we 

m  1  i 

have 


M 


eZ  +  a  Z  -bZ  +cZ 
m,xx  o  m,x  o  m,y  o  m 


-  I 


k-1 


(a.  Z  ,  -  b.  Z  .  +  C.  Z  ,  )  . 

k  m-k,x  k  m-k,y  k  m-k 
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Since  the  initial  conditions  are  Z  =0  on  T,  we  have,  by  induction, 

m 

that  =  0(1)  as  e  ->  0+  uniformly  in  D.  Finally  Z^  satisfies 

(9.5)  in  D  so  that,  with  the  initial  condition  Z..  =  0  on  T  we 

M 

have  that  Z  =  0(1)  as  e  0+  uniformly  in  D.  Thus  we  conclude 
that  the  remainder  term  Z  is  0(1)  uniformly  in  D  as  e  -+  0+. 

4+2N 

As  before,  we  wish  a,  b,  c,  x(s),  y(s)  e  C  so  that 

„  _3+2N,  ,  1+2  (N-n)  _  _  .  XT  j  „  „2 

geC  h  eC  for  n=0,l,...,N  and  Z  e  C  where 

a  . 

n 


10.  An  Extension  of  the  General  Case  of  §5 


As  in  §9,  we  wish  to  have  a,  b,  c,  d  expressed  as  poly¬ 
nomials  of  order  M  in  terms  of  e.  We  shall  construct  a  solution 
of  form 

N  n  N 

u(x,y;e)  =  UQ(x,y)  +  l  e  Un(x,y)  +  u2(x,y;e)  +  u3(x,y;e)  +  e  w(x,y;e). 

n*l 


The  discussion  is  almost  the  same  as  in  §8,  except  solves 

(n  =  1,2, .. . ,N) 


PDE  aU  =  b  U  +CU 
o  n,x  o  n,y  o  n 


-  U 


n-l,xx 


,  ^akUn-k,x 
k«l 


+  e.  U  .  ) 

k  n-k 


in 


-  b.U  , 
k  n-k,y 

(M  +  3M) . 


IC  U  =  0  on  J*. 
n 


* 

For  h  ,  the  initial  condition  becomes  h  =  6  -  U  on  J  . 

n  n  n  n 

Throughout  the  entire  region  D  we  add  to  the  R.H.S.  of  the  equation 
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in  (9.3)  or  (9.4)  the  term 


M 

I  (ak  ^N-hn-k,x  ^N+m--k,y  +  °k  ^N+m-k) 

k=m  * 


(m  =  1,2,...  ,M) 


Note  that  this  does  not  change  the  order  of  Z  (m  =  0 ,1 , . . . ,M-1) . 

m 

Since  the  problem  solved  by  w  in  §8  is  unchanged,  we  have  once  again 
w  =  0(e)  uniformly  in  M  as  e  ->  0+. 


' 


CHAPTER  IV 


J 

WASOW'S  METHOD  FOR  DEALING  WITH  THE  ELLIPTIC  2  ORDER 


PDE  eAu  -  ~  =  f(x,y) 


1 .  The  Statement  of  the  Theorem 

In  [10]  Wasow  proved  the  following  theorem  and  corollaries. 

Theorem.  Let  D  be  a  finite  open  plane  domain  of  finite  connectivity 
bounded  by  a  finite  number  of  arcs  with 
continuously  turning  tangents.  Let  3D 
denote  the  boundary  of  D  and  let  T  be 


*  t 


b 


the  right-hand  boundary  of  D  (see  Fig. 4.1);  and 

T:x  =  q(y)  e  C^([a,b]).  P(x,y)  is  a  point 
P 

in  the  open  domain  D;  PQ(xQ,y)  is  the  point  ^ 

on  T  through  which  the  abscissa  through  P 
passes.  (Note  that  P  may  not  lie  on  the  left- 


X 


Figure  4.1 


hand  boundary  of  D) .  f(x,y)  e  C(D  +  3D)  and  f(x,y)  e  C  (D) .  Consider 


Problem  I 


PDE  |^»f(x,y) 


(x,y)  e  D 


(1.1) 


IC  w(x,y)  =  0 


(x,y)  e  T 


(1.2) 


with  solution 


w(x,y)  =  /  f(C,y)d?  ; 


x 


(1.3) 


Problem  II 


3  u 

PDE  eAu  +  t—  *  f (x,y)  (x,y)  e  D,  e  >  0 

o  X 


(1.4) 


* 
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BC  u(x,y;e)  =  0  (x,y)  e  3D  (1.5) 

with  solution  u(x,y;e). 

Then 

lim  u(x,y;e)  =  w(x,y)  (1.6) 

e  ->  0 

uniformly  in  a  neighborhood  of  the  line  segment  PP^. 

Corollary  1.  If  (1.2)  is  replaced  by  w(x,y)  =  <f>(x,y)  for  (x,y)  £  T, 
where  4>  £  C(T)  and  (1.5)  is  replaced  by 

r  <f>  (x,y)  (x,y)  £  T 

u(x,y ;e)  = 

>•  ^ (x,y)  (x,y)  £  (3D  -  r)  , 

where  is  so  defined  that  u  £  C(3D),  then 

x 

lim  u(x,y;t)  =  <f>  (x  ,y)  +  /  f(£,y)d? 

£  ->  0  x 

o 

uniformly  in  a  neighborhood  of  the  line  segment 
Corollary  2.  If  £  <  0,  then  r  is  replaced  by  3D  -  T. 

2 .  An  Outline  of  the  Technique  that  Wasow  Used  to  Prove  the  Theorem 

2.1.  First  the  proof  of  the  theorem  is  given  for  a  special  case  characterized 
by  the  following  additional  restrictions: 

(A)  Let  the  boundary  3D  of  D  be 

composed  of  two  open  arcs  and 

T  defined  by  the  functions  x  =  ^(y) 
and  x  =  q1(y)  for  a  <_  y  _<  b 
respectively,  and  of  two  closed  line 
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segments  y  =  a  and  y  =  b,  where  it  is  assumed  that 


l0 


for  a  <_  y  <_  b  (see  Fig.  4.2)  . 


(B)  ^(y)  e  C  (^y| a  1  y  1  b})  . 


(C)  f (x,y)  e  C  (D  +  3D) . 
x 


(D)  w(x,y)  =  J  f(£,y)d£  vanishes  on  3D. 


q0(y> 

Define  E(x,y;e)  by  u(x,y;e)  =  w(x,y)  +  e*E(x,y;e).  Then  E(x,y;e) 
solves 

S  E 

PDE  eAE  +  —  -  -  AW  (x,y)  e  D 
3  X 


BC  E  =  0 


(x,y)  e  3D. 


To  show  (1.6),  we  need  only  show  that  E(x,y;e)  is  bounded  in  D  +  3D, 


eAu  + 


If  G(x,y;£ ,n;e)  is  the  Green’s  function  for  the  operator 
over  the  domain  D,  then  E  can  be  expressed  as 


E(x,y;e)  -  /  /  G(x,y ;£ ,n ;e) Aw  d£dn  . 

D 


It  is  shown  that 


G(x,y;£,n;e) 


5-x 

—  K(x,y;C,n;t)  e2e 


(2.11) 


where  K(x,y;£ ,n;e)  is  the  Green's  function  for  the  operator 
2  1 

e  Au  -  u  over  the  domain  D.  By  (i)  the  fact  that 


0  <  K  <  y(r;e)  , 


(2.12) 


where  Y  “  ^  ,  a  Hankel  function  of  the  first  kind  of  order 

zero,  is  a  fundamental  solution  for  Au - -?  u^  ;  and  r  -  /(x-£) ^(y-n) 1 , 

4e 

and  (ii)  the  asymptotic  formulae  of  the  modified  Bessel  and  Hankel 


(1) 


See  [9]  p.  566. 
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functions,  we  can  show  that  E(x,y;e)  is  bounded  in  D. 

2.2.  Removal  of  Assumption  (D) .  Any  domain  satisfying  assumption  (A) 
will  be  called  normal.  Let  D*  by  any  closed  subdomain  of  D  + 

Then  D1  can  be  embedded  in  a  closed  normal  subdomain  D*  of  D  + 
such  that  (i)  the  right-hand  boundary  arc  of  D*  is  part  of 


and  u*(x,y;e)  solves 


PDE 

eAu*  +  =  f* 

3x 

(x,y)  e  D 

BC 

u*  =  0 

(x,y)  e  3D, 

then  by  the  result  in 

§2.1  we  have  that 

lim 

u*  =  pw  uniformly  in 

(D  +  3D)  . 

e  -*■  0 

Let  H(x,y ; e)  =  u(x,y; 

e)  -  u*(x,y;e).  Then 

H  solves 

PDE 

eAH  +  =  f  -  f* 

3x 

(x,y)  e  D 

BC 

H  =  0 

(x, y)  e  3D  . 

X  + 
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To  show  that  (1.6)  in  §1  holds  uniformly  in  D',  we  must  show  that 


lim  H(x,y;e)  =  0  (x,y)  e  D'. 

e  0 


Again 

H(x,y;e)  -  /  /  (f  -  f*)  G(x,y;£,n;e)  d£dn.  (2.21) 

D-D* 


By  (2.11)  and  (2.12)  in  §2.1  and  the  fact  that 


for  —  >  t  ,  where  K  depends  on  the  choice  of 


Since  r  >  r  , 
o 


we  need  only  require  that  e  < 

< 

(£,n)  e  (D  -  D*)  ,  then 


Since  (x,y)  e  D*  and 


r  r  <  r  <  R  , 
° 

0  <  0  <  0. 

L  o  —  —  1  ♦ 


where  R  is  greater  than  the  diameter  of  D.  Thus  -1  <_  cos  0  <  1. 

Let  =  min  |  ,  2tt  -  0^|  (see  Figure  4.4)  , 

(x,y)  e  D' 

(£,n)  e  (D-D*) 


Then 


n  (  .  1  K  "  2e(1  COS  V 

0  <  G(x,y;5,n;e)  <  j  * -  e 


'7T£r 


(2.22)  approaches  to  zero  as  e  ->  0  for  (x,y)  e  D'  and 


(2.22) 
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(£,n)  e  (D  -  D*) .  From  this  and  (2.21),  we  have  that 


lim  H(x,y;e)  =  0  (x,y)  e  D’. 

e  -+  0 

At  this  point  Corollary  2  (i.e.  if  e  <  0,  then  is  replaced  by 

rx)  is  proved.  If  e  <  0,  then  we  replace  y(r;e)  by  y(r;-e). 
Then  asymptotically  we  have 


Y(r;-e) 


y 


so  that 


where  2tt  - 
0  <  a  <  1. 


0  <  G(x,y;£,n;e) 

0^  0  0q,  so  that 


i  K  k  (1+cos  e) 

T  *  - -  e 

/ -7T  e  r 

-a  <  cos  0  <  1  for  some 


a, 


I  00 

2.3.  Removal  of  Assumption  (C) .  Let  {f^ (x,y) } | be  a  sequence 

3 

of  functions  with  f^  e  C  (D  +  3D) ,  which  converges  uniformly  to 
f(x,y)  where  f  e  C^CD)  and  f  e  C(D  +  3D).  Then  we  show  that 

(i)  w^(x,y)  ->  w(x,y)  uniformly  in  (D  +  3D)  as  v  ->  00 , 

(ii)  u^(x,y;e)  ->  u(x,y;e)  uniformly  in  (D  4-  3D)  as  v  and 

(iii)  lim  u  (x,y;e)  =  w  (x,y)  uniformly  in  D*  as  defined  in 

£  +  0  V  V 

§2.2, 


where  w  (x,y)  satisfies 


v 


PDE 


3w 

— ^  =  f 

3x  v 


(x,y)  e  D 


IC 


w  = 

V 


(x,y)  e  T 


0 


and  u  (x,y;e)  solves 


v 


0 
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9u 

PDE 

eAu  +  t  ~  = 
v  9x 

f 

V 

(x,y) 

e  D 

BC 

u  = 

V 

0 

(x,y) 

e  9D. 

From  the  results  (i) ,  (ii) ,  and  (ill)  above,  the  relation  (1.6)  holds 
in  D'  . 


2.4.  Removal  of  Assumption  (B) .  Here  condition  (B)  is  replaced  by 
the  assumption  that  q^(y)  e  Cp  •  Let  be  the  arc  (replacing 


given  by 


such  that 


X  -  q£u)<y>  .  a  1  y  i  b, 


(i)  q^Cy)  e  C3, 

(il)  q0(y)  <  q£u)  <  q0(y)  +  ±  . 


is  the  region  bounded  by  the  arcs  ,  F^ ,  y  =  a  and  y  *  b. 

q£(y)  is  a  function  defined  on  a  <_  y  <_  b 
such  that 


(i)  q2(y)  £  C  , 

(ii)  q1  <  q2  <  qQ,  a  y  ^  b. 

x  -  q2(y) 


Define  a  mapping  x'  =  x  + 


m 


y'  =  y 


where  m  is  a  constant  such  that 


1_ 

y 


0  <  m  <  min  {x  -  q2(y)}  =  min  {qQ(y)  -  q2(y)), 


(x,y)eT 


0 


a<_y<b 


as  shown  in  Fig.  4.5. 
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By  this  mapping  (x,y)  e  (D  +  3D)  is  mapped  into  (x',y')  in  the 
region  bounded  by  the  green  lines  as  x  -  q2(y)  m.  Now  let 


f  (x'  ,y' )  =  f  (x,y) 


f  x'  - 


x'  -  q2(y) 

1  +  mp 


Then  f  e  C^(D^^)  and  f  (x,y)  f(x,y)  uniformly  as  y  -*  00  in  D. 

y  y 

Then  we  show  that 

(i)  w  ■>  w  uniformly  in  D  +  3D  as  y 

y 

(ii)  u^  -*■  u  uniformly  in  D  +  3D  as  y  -*■  00 ,  and 

(iii)  lim  u  =  w  uniformly  in  any  subdivision  of 

e  -*  0  y  y 
x 

where  w^(x,y)  =  /  f^(£,y)d£  solves 


D(y)  +  r(y) 


3w 


PDE  =  f 

3x  y 


(x,y)  e  D 


(y) 


IC 


w  =  0 

y 


(x,y)  e  3D 


(y) 


and  u  (x,y;e)  solves 


3u 

PDE  eAu  +  "r — ~  =  f 
y  3x  ) 


(x,y)  e  D 


(y) 


BC 


u  =  0 

y 


(x,y)  e  3D 


(y) 


By  the  above  results  lim  u(x,y;e)  =  w(x,y)  uniformly  in  D  follows 

e  ->  0 

2.5.  Removal  of  Assumption  (A) ♦  In  this  section  D  is  no  longer 
necessarily  normal.  Let  D’  be  a  subdomain  of  D  +  3D  such  that  D' 
can  be  enclosed  by  a  normal  subdomain  D*  with  its  right  boundary  arc 
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* 

Tq  contained  in  r^.  Let  u(x,y;e)  be  the  solution  of  Problem  II 
in  §1  and  u(x,y;e)  =  $(x,y;e)  for  (x,y)  e  3D*.  In  D*,  u(x,y;e) 
can  be  written  as 

u(x,y;e)  =  u1(x,y;e)  +  u2(x,y;e)  , 
where  u^(x,y;e)  satisfies 

3ul 

PDE  eAu1  +  =  f(x,y)  (x,y)  e  D* 

BC  u^  =  0  (x,y)  e  3D* 

and  u2(x,y;e)  satisfies 

^U2 

PDE  eAu  +  - -  =  0  (x,y)  e  D* 

Z  3  X 

BC  u2  =  $  (x,y)  e  3D*  . 

Since  D*  is  normal  domain,  we  have  by  §2.4  that 

lim  u  =  w 

e  -+  0  i 

uniformly  in  D*.  The  fact  that  (1.6)  holds  uniformly  in  D'  now 

follows  from  the  proof  that  lim  u  =  0  uniformly  in  D’ . 

e  -*■  0 


3 .  Extension 

The  preceding  results  can  be  immediately  extended  to  the 
more  general  partial  differential  equation 


T  r  t  .  ,  3u  ,  3u  _ 

L[u]  =  eAu  +  a  - - 1-  b  —  =  f 

3x  3y 


provided  a  and  b  are  conjugate  harmonic  functions  not  having  a 
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common  zero  in  D  +  3D,  for  then  there  exists  a  conformal  mapping 
taking  L[u]  into 


(a2  +  b2)  {eA'u  +  } 

d  X 


The  above  statement  can  be  shown  as  follows:  Let  w  =  x  +  iy, 

z  =  x'  +  iy'  where  z  =  f (w)  is  an  analytic  function.  Then  we  have 
that  (1)  Cauchy-Riemann  equations  x^  =  y'  ,  Xy  =  ”^x  *  (3.1) 

(2)  Ax’  =  Ay’  =  0.  (3.2) 

Now  by  the  chain-rule  of  differentiation  and  (3.1),  (3.2),  we  have 

eAu  +  au  +  bu  =  e [ (x' )2  +  (x' )2]  (u  ,  ,  +  u  ,  ,) 
x  y  x  y  x  x'  y  y 

+  (ax'  +  bx')  u  ,+  (ay*  +  by')u  i  .  (3.3) 

x  y  x  ^x  y  y 

dz 

If  z  =  f (w)  and  —  =  a  +  ib,  then  a,  b  are  a  pair  of  conjugate 

dw 

harmonic  functions: 


x '  =  a  =  y '  ,  x  ’  =  b  =  -  y '  . 
x  y  y  x 


dz 


We  will  have  a  conformal  mapping  in  D  +  3D  as  long  as  ~  ^0  at 
any  point  in  D  +  3D.  Then  we  have  from  (3.3) 


2  2  2  2 

eAu+au  +bu  =e(a  +b)A’u+(a  +b)ui 
x  y  x 

=  (a2  +  b2)  {eA'u  +  u  ,}  , 

X 


where  A'u=u,  ,+u,  ,  . 

xx  y  y 


CHAPTER  V 


A  NOTE  ON  THE  ASYMPTOTIC  SOLUTION  FOR  HEAT  EQUATION 


In  the  present  chapter  we  apply  the  technique  of  Wasow  [10] 
to  the  heat  equation  in  one  space  variable.  Everything  stated  below 
can  be  carried  over  to  the  case  of  several  space  variables. 


1.  Statement  of  the  Theorem 

Consider  an  open  domain  D  bounded  below  by  the  straight 
line  t  =  0  and  laterally  by  AC  and  CD  where  AC  and  BD  do 
not  include  the  end-points  A  and  B  (see  Figure  5.1).  Consider 


Figure  5.1 


the  partial  differential  equation 


x 


eU  “  U  =  f(x,t)  (x,t)  e  D 

XX  L 


(1. 


with  initial  condition 
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u  =  <|>(x) 

and  the  boundary  conditions'^ 

u  =  ^(xjt)  =  w(x,t) 

u  =  \l>2  (x >  t )  -  w(x,t) 
where  w(x,t)  is  given  by 


e  AB 

(1.2) 

(x,t)  e  AC 

(1.3) 

(x,t)  e  BD 

(1.4) 

w(x,t)  -  4> (x)  -  f(x,x)  dt 


and  satisfies  the  reduced  equation  (i.e.  put  e  =  0  in  (1.1)) 


w  =  -f(x,t)  (x,t)  e  D 

w 


and  the  initial  condition 


w  =  <J>(x)  x  e  AB 


If  f(x,t)  e  C(D+9D)  ,  f(x,t)  e  C^"(D)  and  4>(x)  e  C^AB)  ,  then 


lim  U(x,t;e)  =  w(x,t) 
e  -*•  0+ 


(1.5) 


uniformly 
below  by 


for  (x,t)  e  , 
AB  ,  laterally  by 


where  is  the  closed  domain  bounded 

AC ,  BD  and  above  by  CD  :  t  =  T  . 


(1) 


Note  the  restriction  on  the  boundary  conditions.  This  will  be 
clarified  later. 
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2.  Special  Case 

Now  the  proof  of  the  theorem  will  first  be  given  under  the 
following  additional  restrictions: 

(a)  f(x,t)  e  C2(D+3D)  ,  f(x,t)  e  C3(D) 

(b)  ij)(x)  e  C3(AB) 

Assume 


u(x,y;e)  =  w(x,t)  +  eE(x,t;e) 


(1.6) 


where  E(x,y;e)  solves 


eE  -  E  =  ~w  (x,t)  e  D_ 

xx  t  xx  1 


(1.7) 


with  the  initial  condition 


E  =  0  (x,t)  e  AB 


(1.8) 


and  the  boundary  conditions 


E  =  0  (x,t)  e  AC  and  BD 


(1.9) 


There  exists  a  unique  solution  of  (1.7),  (1.8)  and  (1.9)  of  the  form 


(1) 


f  r 

E(x,t;e)  =  |  |  G(x,t;£,i;e)w  d?  dx 

JDi J  “ 


(1.10) 


where  G  = 


G  -  v  and 
o 


G  (x,t;C,i;e) 
o 


2/tte  (t-x) 


(x-o 

4e (t-x) 


(1) 


The  proof  can  be  found  in  [9],  pp.  521-525;  209-210. 
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is  the  fundamental  solution  of  the  heat  equation 

eU  -  U  =  0  (x,t)  e  D-  . 

xx  t  1 

The  term  v(£,T;x,t;e)  ,  when  considered  as  a  function  of  £  and  t 
containing  the  parameters  x,t,e,  satisfies 

eV_r  +  V  =0  for  t  <  t  (x,t)  e  D. 

££  t  1 

with  the  initial  condition 

V  =  0  on  CD  :  t  =  t 

V  =  G  (x,t;£,r;e)  on  AC  and  BD 

o 

and  is  of  class  C^  in  the  interior  of  D  ;  or  v(x ,  t ;  £  ,t  ;e  )  satisfies, 
with  resnect  to  x  and  t  , 

eV  ~  V.  =  0  for  t  >  t 

xx  t 

and  V  =  0  for  t  =  t  (1.11) 

V  =  G  on  AC  and  BD  , 

o 

V  and  G  denend  only  on  the  difference  t  -  x  .  Note  that  as  t  ■+  t  +  0  , 
with  fixed  x  not  equal  to  £  ,  the  boundary  values  (1.11)  are  continuous 
at  the  corners  of  tne  region.  On  observing  that  the  boundary  values 
are  non-negative,  we  have  by  the  maximum  principle  for  the  parabolic 
equation^)  that  V  >  0  and  G  _>  0  in  D^  .  Hence  by  the  fact  that 


(1) 


See  [A],  Chapter  II. 
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G 


=  G 


o 


V 


> 


G(x,t;£,x;e)  _<  Go(x,t;£,x;e) 


in  D1+9D1 


Now  we  shall  show  that  E(x,t;£,x;e)  is  bounded  in  . 

1  2 

By  the  restrictions  (a)  and  (b),  we  have  w  e  C  (D.. )  and  E  e  C  (D, )  . 

XX  X  X 

By  the  fact  that  0  £  G  _<  Gq  it  suffices  to  find  an  upper  bound  for 


We  recall 


f 

J  J 

D, 


G  (x,t;£,x;e)d£  dx 
o 


GQ(x,t;5,T ;e) 


(x-er 

Ac (t-x) 


2/rre(t-x) 


Thus 


G  dg  dx  = 

•  /• 

■*  D,  • 

°  J 

D1 

_  (x-gT 

Ae (t-T )  ,  , 

e  dg  dr 


T  c°° 


o  -1  -°°  2/jTe(t-x) 


_  (x-O 

e  4E(t"T)  d5  dT 


T  r°°  1  2 

1  “CL  j  j 

-  e  da  dx 


J  o  ^  -  0°  /rr" 


where  we  have  let 


x~g 


2 v  e(t-x) 


w  Si  i — 

=  a  .  Then  by  the  fact  that  /  e  da  =  vtt  , 


we  obtain 
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r  r 

I  G  d£  dx  < 

J°1  J  ° 


dt  =  T 


(1.12) 


From  (1.10)  and  (1.12),  we  get 


E(x,t:e)|  = 


= 

1 

f  <• 

J 

D-  * 

1 

G(x,t;£,x;e)w  ^d?  dx 


Di J 


G  (x,t;C»T;e)  w  d£  dx 

o 


<  MT 


(1.13) 


where  M  =  max  |w  J  .  Hence  it  follows  from  (1.6)  and  (1.13)  that 
in  D  XX 

lim  u(x,t;e)  =  w(x,t) 
e  -*•  0 

uniformly  for  (x,t)  e  . 

Note  that  if  in  (1.3)  and  (1.4)  t  w  and  ^  w  on 

AC  and  BD  respectively,  i.e.  arbitrary  prescribed  boundary  conditions, 
then  (1.10)  will  be  replaced  by 


E(x,t;e)  = 


AC 


n  \  3G  a 
«rw)  dT  - 


BD 


(*2'w)  If  dl 


G  d£  dx  . 


(1.14) 


Then  in  order  to  find  the  upper  bound  of  E(x,t;e)  ,  we  not  only  have 


to  check  the  last  integral  of  (1.14)  as  we  have  done  above,  but  also 

3G  3g  8Go 

f  —  dx  and  /DT.  —  dx  .  By  the  same  token  we  check  f.  — rr-  dx 

AC  dt,  BD  dt,  AL  dt, 


dx  instead. 


and 


3G 

o 


dx 


m  _  C  X~g  ) 

x-g _  e  4e  (t-i) 

Jo  4/rT[e(t-x) 


dx 


(1.15) 


Let 


x-e 


2/sCt-x ) 


=  3 


Then  d3  = 


x-g 


4*/e^(t-x) 


3/2 


dx  and  (1.15)  becomes 


x-g 

,  2/e  (t-T)  _ft2 

—  e  P  dS  (1.16) 

e/rr"  ^  x-g 

2 /it” 


By  the  fact  that 


lim 
e  ->  o 


x-g 

•  2/ e (t-T) 

x-g 

2/e"t 


This  with  (1.15)  and  (1.16) ,  we  get 


— r^r  dx  =  0(— )  as  e  ->  0 

Uc  35 


Similarly 


3G 

— ~r  dx  =  0  (  )  as  e  ■>  0 
J  BD  3£ 

In  this  case  E  is  not  bounded  and  we  get  a  boundary  layer  phenomenon 
along  AC  and  BD  ,  which  are  the  characteristics  of  the  reduced 


equation.  S.L.  Kamenomostskaya  [6]  shows  that  on  AC  and  BD  there 
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is  a  boundary  layer  term  z  such  that  |z|  <_  K  e 


-(g.//e) 


+  /e  Z  , 


>  0  (i  =  1,2)  near  AC  and  BD  ,  and  |z|  K  e  +  /e  Z 


6  >  0  ,  elsewhere.  Also  z  =  u  -  w  on  AC  and  BD  . 


(1) 


3.  Removal  of  Restrictions 


Now  we  shall  relax  the  restrictions  (a)  and  (b)  in  §  2  to 

(a) ’  f  e  C(EH-3D),  f(x,t)  e  C1(D) 

(b )  *  <f>(x)  £  C^(AB) 


as  in  the  statement  of  the  theorem. 


Let  {f  (x,t)}  ,  be  a  sequence  of  function  with  f  e  C 

V  'v=l  n  V 

3 

and  f  e  C  (D)  such  that  f  ->  f  uniformly  in  D+3D  as  v  -*  1:0  . 

v  v 

i 00  3 

Let  {d>  (x)}  ..  be  a  sequence  of  functions  with  <t>  eC  (AB)  such 

V  1 V=1  '  V 

that  <f>  ->  <J>  uniformly  in  AB  as  v  ->  00  .  Let  w  (x,t)  solve 


PDE 


w  =  -  f 
t  v 


in  D^+3D^ 


IC 


w  =  <f> 


on  AB 


Then 


w  =  4> 


v  v 


f  (x,T)dr 


V 


Let  u  (x,t;e)  solve 


v 


(1) 


This  result  is  similar  to  V  in  §6  ,  Chapter  III. 


(D+3D) 
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PDE 

eAu  -  u  = 
t 

f 

V 

in  D^+3D^ 

IC 

u  =  d> 

V 

on  AB 

BC 

e 

ii 

K 

c 

on  AC  and 

Then  (i) 

w  ->  w 

V 

uniformly 

in 

(D^+3D^)  as  v 

(ii) 

u  ->  u 

V 

uniformly 

in 

(D^+3D^)  as  v 

(iii) 

lim  u 

e  -*  o 

^  uniformly  in  . 

Result  (i) : 


w  -  w 

V 


i  If,  - 


(f  -f)dt 


and 


<  A  +6  •  T  ,  for  v  >  N 

—  v  v 


where  A 

v 


max  |  <f>  —  4> 

(x,t)  e  AB 


and  6  = 

v 


max 

(X,t)  e  D1+8D1 


f 


lim  w 
v  ->  00 


w 

v 


uniformly  in  D^+9D^ 


Result  (ii) : 


Let 

V  =  u  -  u  . 

V 

Then 

V  solves 

PDE 

eV  -  V 

-  f 

-  f  in  D.+8D 

XX  t 

V 

1 

IC 

i 

-e-" 

II 

> 

on  AB 

BC  V  =  w  -  w 

v 


on  AC  and  BD  . 
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Then 


V  = 


J  AB 


(<4>  -  d>)G(x,t  ;£ , t  ;  e)d£  +  (w  -  w)  77  dx 


AC 


v  y  9? 


+ 


BD 


(w  -  w)  ~  dt  + 

V  0  c, 


JD,  j 


(f  -  f ) • G  d£  di 


v 


For  fixed  e  >  0  ,  we  have  that 


lim  V  =  0 

v  -*  00 


For  this  particular  v  result  (iii)  follows  from  the  nroof  in  §2  . 

Now  to  show  that  u  -*  w  uniformly  as  e  -*  0  for  (x,t)  e  , 

we  must  show  that  for  every  6  >  0  ,  there  exists  a  e  =  e(6)  such 
that  0  <  e  <  e  implies  that 

|u  -  w  |  <  6  (x,  t)  e  D 

So  let  us  choose  6  >  0  .  Then 

(i)  choose  N^(6)  sufficiently  large  that  v  >  N^  implies 
that 

|w^  -  wj  <  6/3  (x,t)  e  , 

(ii)  choose  N£(6)  sufficiently  large  that  v  >  N^  implies 
that 


1  ’ 


u 


<  6/3 


(x,t)  e  D 


-  84 


(iii)  choose  e(6)  such  that  0  <  e  <  e  implies  that 


u  -  w  <  6/3  (x, t)  e  D- 

1  V  V  1  1 

Hence  for  v  >  N  =  maxiN^jN^}  and  0  <  e  <  e  we  have 


u  -  w  <  u  -  u  +  u  -  w  +  w  -  w  <  6 

'  —  1  v  1  1  v  v1  1  v  1 


(x,t)  £  D1 


CHAPTER  VI 


SOME  QUESTIONS  AND  PROBLEMS 

In  conclusion  we  formulate  some  questions  and  problems  which 
naturally  arise  from  the  preceeding  considerations. 

(1)  We  assume  the  functions  in  hypothesis  in  Chapter  II 

6  4 

are  all  C  and  in  hypothesis  in  Chapter  III  are  all  C  .  We 

could  use  Wasow's  technique  which  is  presented  in  Chapter  IV  to  relax 

these  restrictions. 

(2)  The  problem  of  extending  these  methods  to  the  same  types 

th 

of  partial  differential  equations  of  n —  order  in  n  dimensions  for 
n  >  2 . 

(3)  It  is  desirable  to  pursue  an  investigation  of  the  asymptotic 

solution  u  of 

e 

L  [u]  =  e  (u  -  u  )  +  au  +  bu  +  cu  =  d 
e  xx  yy  x  y 

in  R,  a  region  bounded  by  two  characteristics  of  L£[u]  =  f  and  two 
characteristics  of  the  reduced  equation 

L  [u]  =  aU  +  bU  +  cU  =  d  , 
o  x  y 

as  shown  in  Figure  6.1  below,  with  the  prescribed  values  u^  =  U  on 

I  and  u  ^  U  on  J.  It  is  clear  that  we  should  have  a  quite  different 
e 

maximum  principle  for  L  [u]  =  f  from  the  maximum  principles  that  we 


. 
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have  for  the  elliptic  and  parabolic  equations. 

Then  there  naturally  arises  the  question:  "What  happens 
when  the  characteristics  of  higher  order  equation  and  the  characteristics 
of  the  lower  order  equation  coincide?" 
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